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Abstract: With the goal of identifying signatures that select specific neutrino 
oscillation parameters, we test the robustness of global oscillation solutions that fit 
all the available solar and reactor experimental data. We use three global analysis 
strategies previously applied by different authors and also determine the sensitivity 
of the oscillation solutions to the critical nuclear fusion cross section, 5*17 (0), for the 
production of *B. Our standard results make use of the precise new measurement of 
5*17(0) by Junghans et al. The globally favored solutions are, in order of goodness of 
fit: LMA(the only solution at 2cr), LOW, and VAC. The NC to CC ratio for SNO 
is predicted by the standard global analysis to be 3A5~^Ql^{la) which is separated 
from the no-oscillation value of 1.0 by much more than the expected experimental 
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error. The predicted range of the day-night difference in CC rates is 8.31^:^(1(7 
A measurement by SNO of either a NC to CC ratio > 3.3 or a day-night difference 
> 10%, would favor a small region of the currently allowed LMA neutrino parame- 
ter space. The global oscillation solution predicts a ^Be neutrino-electron scattering 
rate in BOREXINO and KamLAND in the range 0.65l|]:[J^(la) of the BPOO standard 
solar model rate, a prediction which can be used to test both the solar model and 
the neutrino oscillation theory. Only the LOW solution predicts a large day-night 
effect(< 42%, 3a) in BOREXINO and KamLAND. For the reactor KamLAND ex- 
periment, the LMA solution predicts a charged current rate relative to the standard 
model of 0.441q;q7(1(t), -^threshold = 1-22 MeV. We have also evaluated the effects of 
including preliminary Super-Kamiokande data for 1496 days of observations. 
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1. Introduction 

How robust are the predictions of the allowed neutrino oscillation solutions that 
are obtained by fitting the currently available solar neutrino data with different 
neutrino oscillation parameters? How can we determine experimentally the neutrino 
oscillation parameters? We address these questions in the present paper. 

To determine the robustness of the oscillation solutions and their predicted ex- 
perimental consequences, we compute the currently allowed regions in neutrino os- 
cillation space using three different analysis strategies, each strategy previously ad- 
vocated by a different set of authors. We include all the available solar neutrino 
and reactor oscillation data [|l], H, H, ^, |, |], ^ . We also present "Before and After" 
comparisons of the globally allowed solutions using both the previous(1998) standard 
value of 5*17(0), the low energy cross section factor for the production of ^B, and a 
precise new measurement for this important quantity. 

In order to identify which quantities are most likely to lead to an experimental 
determination of the neutrino oscillation parameters, we use the allowed regions in 
neutrino parameter space to predict the expected range of the most promising quanti- 
ties that can be measured accurately in solar neutrino experiments. For the Sudbury 
Neutrino Observatory (SNO) |[, we report the predictions of the currently allowed 
regions for the neutral current to charged current double ratio, [NC]/[CC], the day- 
night effect, and the first and second moments of the electron recoil energy spectrum. 
For experiments like BOREXINO [l^ and KamLAND that can detect ^Be solar 



neutrinos, we evaluate the predicted range of the neutrino-electron scattering rate 
and the difference between the night and the day event rates. 



It pays to be lucky. We show in figure |TT] and figure |l^ that if SNO measures (as 



predicted by some LMA oscillation parameters) a large value for either the neutral 
current to charged current ratio or the day-night difference in rates, or if BOREXINO 
or KamLAND measures(as predicted by some LOW oscillation parameters) a large 
day-night difference, then this measurement will uniquely select the correct neutrino 
oscillation solution and fix the neutrino oscillation parameters to within rather small 
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uncertainties. On the other hand, a small value of the neutral current to charged 
current ratio would favor the SMA solution. 

How should this paper be read? Since you got this far, we recommend that you 
now turn to the final section, the summary and discussion section, and see what 
results seem most interesting to you. The summary and discussion section refers to 
what we believe are the most useful figures and tables, which the reader may wish 
to glance at before deciding whether to invest the time to read the detailed text in 
the main sections. 

This paper is organized as follows. In section ^ we present the predicted so- 
lar neutrino fluxes and their uncertainties that are derived from the BPOO standard 



solar model |T2[ with and without including the recent precise measured value for 
5*17(0) determined by Junghans et al. p. We determine in section |^ the regions in 
neutrino oscillation parameter space that are allowed using three different analysis 
scenarios; we present results obtained using both the new value for 5*17(0) and the 
previous(1998) standard value. We use the three sets of global oscillation solutions 
in section ^ to determine for SNO the 3a (and la) ranges of predicted values for the 
neutral current to charged current ratio, the difference between the day and the night 
CC event rates, and the first two moments of the energy spectrum of the CC recoil 
electrons. In section ^, we determine for the BOREXINO and KamLAND solar neu- 
trino experiments the 3a (and la) ranges of allowed values for the neutrino-electron 
scattering rate and the day-night effect. For the reactor KamLAND experiment, we 
present in section ^ the predicted 3a (and la) ranges of the charged current (an- 
tineutrino absorption) rate and the distortion of the visible energy spectrum. We 
summarize and discuss our conclusions in section |^ 

We use, unless stated otherwise, the techniques and parameters for the analysis 
that we have described elsewhere O, |I4|, O, llH]. 



2. BPOO + New 

In this section, we present the predicted solar neutrino fluxes (table |1|) and their 
uncertainties (table 0) that were computed with the BPOO solar model |T^, with 
and without taking account of a recent precise measurement of the production cross 
section that leads to the emission of ^B neutrinos. Both the solar neutrino fluxes 
and their uncertainties are used in calculations of the allowed regions for neutrino 
oscillation parameters. 

There has been a lot of beautiful experimental work devoted to measuring the 
rate of the fusion reaction ^Be(p, 7)^B p, 17|. Recent experiments have concentrated 



on reducing the systematic uncertainties, in going to lower center-of-mass energies 
that are more relevant to the sun, and in using new measurement techniques. The 
goal is to reduce the uncertainty in the low energy cross section factor for the reaction 
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Be(p, 7) B to less than or equal to 5% |18|, so that it is no longer the dominant un- 



certainty in the prediction of the solar neutrino flux. Very important experiments 
are continuing to be performed on this reaction. 

The recent Junghans et al. precision measurement of the low energy cross section 
for the fusion reaction ''Be(p,7)^B yields 

Sn{0) = (22.3 ± 0.7 (expt) ± 0.5 (theor)) eVb. (2.1) 

The quoted uncertainty for the Junghans et al. measurement, which combines sta- 
tistical and systematic uncertainties, is smaller than for all previous low energy mea- 



surements. Other relatively recent measurements [0] have yielded values of 6*17(0) 
that are smaller than the Junghans et al. value. 

For illustrative purposes, we adopt in this paper 5*17(0) = (22.3 ± 0.9) eVb in- 
stead of the previous (1998) standard value of 5*17(0) = 19l26Vb [|19[ that was used 
in the original calculations of the BPOO standard model neutrino fluxes. The 
value of 5*17(0) given in eq. ( |2.1| ) is 0.8cr larger than the previously used value. 

We will present throughout this paper the differences in the calculated quantities 
that result from using ^17(0) = (22.3 ± 0.9) eVb instead of ^17(0) = 19+^eVb. These 
two values of 5*17(0) span the set of recent experimental results pjf . 

2.1 Fluxes 

Any plausible change in the value of 5*17(0) has no discernible numerical effect on 
the standard solar model |18|. The reaction ''Be(p,7)'^B is so rare (it occurs of order 
twice in every 10^ completions of the p — p chain) that even a factor of ten change in 
the cross section does not affect, to the required accuracy, the calculated solar model 
temperatures and densities. The 15% increase in the best-estimate value of 5*17(0) 
represented by eq. ( |2.1| ) does not change signiflcantly any of the computed physical 
variables in the standard solar model except the ^B neutrino flux. 

For a given solar model with specifled distributions of temperature, density, and 
composition, the calculated ®B neutrino flux depends linearly upon the adopted 
value of ^17(0), i.e., (/.(^B) oc 5i7(0). 

Table |l| presents the predicted neutrino fluxes, and the predicted chlorine and 
gallium capture rates, for the BPOO standard solar model and the ^Be(p,7)®B cross 
section factor of eq. ( p.lj) . Except for 0(^B), all of the fluxes in table are the same 
as in table 7 of the original BPOO paper |jT2|. 

2.2 Uncertainties 

Table ^ summarizes the uncertainties in the most important solar neutrino fluxes 
and in the CI and Ga event rates that are due to different nuclear fusion reactions 
(the flrst four entries), the heavy element to hydrogen mass ratio (Z/X), the radia- 
tive opacity, the solar luminosity, the assumed solar age, and the helium and heavy 
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Source Flux CI Ga 

pp 5.95 X 10° (1.00+°;°}) 

pep 1.40 X 10-2 (i.ooi[5:[!}^) 

hep 9.3 X 10-^ 

^Be 4.77 X 10-^ (1.001°-}°) 

SB 5.93 X 10-4 (1.001°-}^) 

i^N 5.48 X 10-2 (l.00l°:2i) 
4.80 X 10-2 (l.00l°-2^) 

i^F 5.63 X 10-4 (l.00l°:25) 



(SNU) 


(SNU) 


n nn 

U .UU 


fiQ 7 


0.22 


2.8 


0.04 


0.1 


1.15 


34.2 


6.76 


14.2 


0.09 


3.4 


0.33 


5.5 


0.00 


0.1 



Total 8.59+};^ 130+^ 

Table 1: Standard Model Predictions (BPOO + New ^B): solar neutrino fluxes and 
neutrino capture rates, with la uncertainties from ah sources (combined quadratically) . 
The neutrino fluxes are the same as in the orig inal BPOO model @ except for the ^B 
flux, which is increased because of the larger adopted value of 5i7(0), see eq. (^?l|). Using 
the 1998 standard value £'17(0) = 19l^eVb 0, the ^B neutrino flux was calculated previ- 
ously to be (f) (^B) = 5.05 (l.00l°;2°). The total rates were calculated using the neutrino 



absorption cross sections and their uncertainties that are given in ref. |20| 



element diffusion coefficients. In addition, the + p reaction causes a 0.2% uncer- 
tainty in the predicted pp flux, a 0.1 SNU uncertainty in the CI event rate, and a 1 
SNU uncertainty in the Ga event rate. 

The uncertainty in the laboratory measurement of the low energy cross section 
factor for the ^He(4He,7)^Be reaction is the dominant uncertainty in the calculation 
of the ''Be neutrino flux and one of the two largest uncertainties in the calculation 
of the ^B neutrino flux. 



For the original BPOO solar model predictions |]T2|, the uncertainties in the in- 
dividual parameters and in the fluxes are, with the exception of the two non-zero 
entries in the column under ^Be -|- p, the same as in table ^ Instead of the entries 
0.040 under ^Be + p, the original BPOO model predictions had a much larger uncer- 
tainty, 0.105. The original BPOO uncertainties were used, for example, in the analyses 

o, mm 



of solar neutrino oscillations described in refs. 

The input data that were used to calculate the uncertainties listed in table 
are given in the papers in which the BPOO model |1^ and the BP98 model [23 
were originally presented. The principal change from the uncertainties that were 
calculated for the BP98 model (cf. table 2 in ref. |^), and which were used in 
many previous analyses of solar neutrino oscillations, is that, assuming the validity 
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< Fractional 


DP 


^He^He 


^He^He 


'^Be + p 


z/x 


opac 


lum 


age 


diffuse 


uncertainty > 


0.017 


0.060 


0.094 


0.040 


0.061 




0.004 


0.004 


0.15 


Flux 




















PP 


0.002 


0.002 


0.005 


0.000 


0.004 


0.003 


0.003 


0.0 


0.003 


^Be 


0.016 


0.023 


0.080 


0.000 


0.034 


0.028 


0.014 


0.003 


0.018 


8B 


0.040 


0.021 


0.075 


0.040 


0.079 


0.052 


0.028 


0.006 


0.040 


SNUs 




















CI 


0.3 


0.2 


0.6 


0.3 


0.6 


0.4 


0.2 


0.04 


0.3 


Ga 


1.3 


1.0 


3.3 


0.6 


3.1 


1.8 


1.3 


0.20 


1.5 



Table 2: For the standard solar model (BPOO + New *B), the average uncer- 
tainties in neutrino fluxes and event rates due to different input data. The flux 
uncertainties are expressed in fractions of the total flux and the event rate uncertainties 
are expressed in SNU. The ''Be electron capture rate causes an uncertainty of ±2% [24| 
that affects only the ^Be neutrino flux. The average fractional uncertainties for individual 
parameters are shown. 



of eq. (|2.1| ), the uncertainty in the cross section for the ''Be(p, 7)*B reaction is reduced 
by a factor of more than three and the estimated uncertainty in the heavy element to 
hydrogen ratio, Z/X, is increased by almost a factor of two( cf. discussion in section 
5.1.2 of ref. [0). 

The predicted event rates for the chlorine and gallium experiments make use of 
neutrino absorption cross sections from refs. [EDI. The uncertainty in the prediction 



for the standard (no oscillation) gallium rate is dominated by uncertainties in the 
neutrino absorption cross sections, 7 SNU (5% of the predicted rate). The uncertain- 
ties in the chlorine absorption cross sections cause an error, ±0.22 SNU (2.6% of the 
predicted rate), that is relatively small compared to other uncertainties in predicting 
the rate for this experiment. For calculations that involve non-standard neutrino 
energy spectra that result from neutrino oscillations or other new neutrino physics, 
the uncertainties in the predictions for currently favored solutions (which reduce the 
contributions from the least well-determined ^B neutrinos) will in general be less 
than the values quoted here for standard spectra and must be calculated using the 
appropriate cross section uncertainty for each neutrino energy 



3. Global neutrino oscillation solutions with different analysis 
prescriptions 

We present in this section the global solutions for the allowed neutrino oscillation pa- 
rameters that are derived using the neutrino fluxes and uncertainties given in table |I| 



- 6 - 



and table We describe in section |3.1| the oscillation solutions that are allowed with 
three different analysis prescriptions (see figure [^) and show how adopting the new 
Junghans et al. value of S'i7(0) reduces the allowed regions in neutrino parameter 
space (see figure ^. We have at different times used all three of the analysis strate- 
gies and various colleagues have advocated strongly one or the other of the strategies 
described here. In section pj.3|, we describe in more detail the three different analysis 



prescriptions. Section Wl. is intended for a general readership, while section is 
intended primarily for aficionados of neutrino oscillations. 

The interested reader may wish to consult also a number of recent papers, 



refs. [[161 , |21| , |2^ , P5[ |26| , |27| , that have determined the allowed solar neutrino os- 
cillation solutions including the CC data from SNO. We will spare the reader erudite 
comparisons, in the cases where there is overlap in the calculated quantities, between 
our detailed results and those of the other authors referred to above. In general, when 
the same analysis procedures and input data are used, all authors obtain consistent 
results. 

In this paper, we study the sensitivity of the calculated neutrino predictions 
to the adopted analysis strategy by comparing the numerical results obtained with 
three different analysis procedures. We also compare, for global neutrino solutions 
and for predicted neutrino measurable quantities, the results that are obtained using 
the previous standard value for S'i7(0) with the results that are obtained using the 
recent precise measurement by Junghans et al. of S'i7(0). In our previous paper on 
solar neutrino oscillations [0, we used only the analysis strategy (a) (cf. Section 3.3 
for a description of all three analysis strategies) and the previously standard value 
for 5*17(0). We apply the sensitivity studies described in section 3 to the calculation 
of a variety of solar neutrino measurables in sections 4-6, providing for the first 
time a quantitative evaluation of the robustness of the predictions to diverse analysis 
strategies. 

3.1 Global oscillation solution: three different analyses 

We first discuss in section |3.1.1| the new global solutions for each of three different 
analysis strategies and then compare in section |3.1.2| the global solutions that are 



obtained with the older and newer (cf. eq. |2?T|) values of 5*17(0). 
3.1.1 New global solution 

Figure shows the allowed ranges of the neutrino oscillation parameters, Am^ and 
tan^ 6, that were computed following three different analysis approaches that have 
been used previously in the literature. These approaches are described in more detail 
in section |3.3| . We follow refs. [28, 29[ in using tan^ 6' (rather than sin^ 2^) in order 
to display conveniently solutions on both sides of 9 = tt / 4. 

Figure |l]a presents the result for our standard analysis, which includes the BPOO 
predicted fluxes and uncertainties and all the experimental data except the Super- 
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Figure 1: Global neutrino oscillation solutions for three different analysis 
strategies. The three different analysis strategies are described in section 3.3. In all 
cases, the value of 5'i7(0) reported in ref. and shown in eq. was used. The input 
data include the neutrino fluxes and uncertainties predicted by the BPOO solar model |jl2| 
and the total measured event rates from the SNO Q, Chlorine and Gallium (aver- 
aged) ^ m, 1^ experiments, as well as the recoil electron energy spectrum measured by 
Super-Kamiokande Q during the day and separately the energy spectrum measured at 
night. The C.L. contours shown in the figure are 90%, 95%, 99%, and 99.73% (3cr). The 
allowed regions are limited by the Chooz reactor measurements to mass values below ~ 
(7-8) xlO~^eV^. The local best-fit points are marked by dark circles. 



Kamiokande total event rate. The measured day and night electron recoil energy 
spectra together represent the Super-Kamiokande total rate in this approach. Fig- 
ure 0b was constructed using an analysis similar to that used to construct figure |l|a 
except that for figure |I]b the total Super-Kamiokande rate is included explicitly to- 
gether with a free normalization factor for the day and night recoil energy spectra. 
Figure |l]c shows the results for the free analysis, in which the calculation is the 
same as for the standard case, figure |I|a, except that the neutrino fiux is not 
constrained by the solar model predictions. 

The allowed regions obtained in the three panels of figure although different, 
do not depend dramatically on the adopted strategy. The allowed regions are reduced 
in size for strategy (a) and strategy (b) as compared to the regions derived using the 
1998 standard value of 5*17(0) |]19|, cf. the results presented in refs. [|T6] and [p^] , 
respectively. 

The most restricted allowed regions are obtained for strategy (b), which has been 
used by the Bari group As discussed in ref. [0, the smallness of the allowed 
regions for this strategy results from the absence in the analysis of the ^B theoretical 
error for the Super-Kamiokande spectrum normalization. We will discuss this point 
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more in the next section. 

The largest allowed regions in figure |^ correspond to strategy (c), used among 
others by Krastev and Smirnov in which the neutrino fiux is treated as a free 
parameter. This was not the case when using the 1998 standard value of S'i7(0 |[T9[| ) 



(see fig. 9 in ref [|1^]). The only changes in the quantities derived from BPOO and 
BPOO + New '^B are the ^B neutrino fiux and its associated uncertainty. Therefore 
the allowed regions are unchanged when going from BPOO to BPOO + New ^B in 
strategy (c), for which the ^B neutrino fiux is treated as a free parameter. Thus 
Fig. [^c is equal to the right panel of Fig. 9 in Ref. |]16| and it is also in excellent 
agreement with the results from the free ^B analysis of ref. |2l|. The change in 
the relative size of the allowed regions between strategy (a) and (c) is due to the 
reduction in size, discussed above, of the allowed regions for strategy (a). 

In constructing figure |I], we assumed that only active neutrinos exist. We derive 
therefore the allowed regions in using only two free parameters: Am? and tan^ 9. 
As discussed below, oscillations into pure sterile neutrinos are highly disfavored (see 
table H) with the larger value of 6*17(0) given in eq. (p.l[), and are not included in our 
standard demarcation of the allowed regions. The allowed regions for a given CL. 
are defined in this paper as the set of points satisfying the condition 

x'(Am^ 6) - xLn < Ax'(C.L., 2 d.o.f.), (3.1) 

with Ax^(C.L., 2 d.o.f.) = 4.61, 5.99, 9.21, and 11.83 for CL. = 90%, 95%, 99% 
and 99.73% (3cr). We use the standard least-square analysis approximation for the 
definition of the allowed regions with a given confidence level. As shown in ref. 
the allowed regions obtained in this way are very similar to those obtained by a 
bayesian analysis. 

Table || gives for our standard analysis strategy (cf. figure |l|a) the best-fit values 
for Am^ and tan^ 9 for all the neutrino oscillation solutions that were discussed in 
our previous analysis in ref. |jT6|. The table also lists the values of Xmin each 



solution. The regions for which the local value of Xmin exceeds the global minimum 
by more than 11.83 are not allowed at 3(T CL. Thus, all solutions with Xmm bad or 
worse than the VAC Sterile solution are not allowed at 3cr for the standard analysis. 

The set of allowed solutions at 3cr C L. would not change if we considered active- 
sterile admixtures and defined the allowed regions as in figure 1 of ref. [|1^] in terms 
of Ax^(CL., 3 d.o.f.) instead of Ax^(CL., 2 d.o.f.). Table ^ shows that all regions 
with Ax^ > 11.83 (corresponding to 3(7 CL. for 2 d.o.f.) also satisfy Ax^ > 14.16 
(corresponding to 3(T CL. for 3 d.o.f.). 

Many authors(see, e. g., ref. and references quoted therein) have discussed 
the possibility of bi-maximal neutrino oscillations, which in the present context im- 
plies tan^ 6* = 1. Figure |l| shows that bi-maximal mixing is disfavored for the pre- 
ferred LMA solution. Quantitatively, we find that there there are no solutions with 
tan^ = 1 at the 99.87% CL. for our standard analysis strategy, corresponding to 
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Solution 






tan^(^) 




2 

Xmin 


g.o.f. 


LMA 


3.7 


X 


10- 


-5 


3.7 X 10- 


-1 


34.5 


67% 


LOW 


1.0 


X 


10- 


-7 


6.7 X 10- 


-1 


40.8 


39% 


VAC 


4.6 


X 


10- 


10 


2.5 X 10 





42.3 


33% 


Sterile VAC 


4.7 


X 


10- 


10 


3.0 X 10 





49.1 


13% 


SMA 


5.2 


X 


10- 


-6 


1.8 X 10- 


-3 


49.9 


11% 


Just So^ 


5.5 


X 


10- 


12 


0.61(1.6) X 


10° 


52.1 


7.8% 


Sterile Just So^ 


5.5 


X 


10- 


12 


0.61(1.6) X 


10° 


52.1 


7.8% 


Sterile SMA 


4.6 


X 


10- 


-6 


3.4 X 10- 


-4 


52.3 


7.5% 



Table 3: Best-fit global oscillation parameters with all solar neutrino data. 

The table gives for the global solution illustrated in figure |^a the best-fit values for Am^, 
tan^0, Xmin' ^^'^ g-O-f- foi" s-ll the neutrino oscillation solutions that were discussed in our 
previous analysis in ref. [^]. The differences of the squared masses are given in eV'^. The 
number of degrees of freedom is 39 [38(19 night spectrum, 19 day spectrum) + 3(rates) 
— 2(parameters: Am^, 9)]. The BPOO best-fit fluxes and their estimated errors have been 
included in the analysis (see table |l] and table ^). The rates from the GALLEX/GNO and 
SAGE experiments have been averaged to provide a unique data point. The goodness-of-flt 
given in the last column is calculated from the value of x^/d.o.f at each local minimum 
(i.e., for LMA, SMA, VAC). Solutions that have Xmin ^ 34.5 + 11.83 are not allowed at the 
3cj C.L. 



panel (a) of figure |1]. For the strategies corresponding to the panels (b) and (c), 
respectively, there are no solutions with tan^ 6^ = 1 at the 99.95% C.L. and 98.95% 
confidence level. Maximal mixing is allowed for the less favored LOW and VAC 
solutions at a C.L. that varies between 95% and 3a, depending upon the analysis 
strategy. In summary, maximal mixing is not favored but it is not rigorously ruled 
out. 

The upper limit on on the allowed value of Am^ is important for neutrino oscil- 
lation experiments, as stressed in ref. |21]. In units of 10~^ eV^. we find the following 
3a upper limits on Am^ (values in parenthesis are obtained with the older value of 
5*17(0), see ref. Il^): case a: 4.7 (7.5); case b: 3.0 (6.8); and case c: 6.5 (6.5). The 



upper limit therefore lies between 3.0 and 7.5, in units of 10 ^ eV^, depending upon 
what is assumed about 5*17(0) and the analysis strategy. 



3.1.2 Global "Before and After" 

Figure ^ compares the allowed regions found with the new theoretical neutrino 
flux and uncertainties (table and table ^ with our previously published results(see 
figure 9a of , ref. |]16|)- We like to refer to this comparison as our global "Before- After" 
figure. 
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Figure 2: Global "Before and After." The left panel shows the allowed regions for 
the neutrino oscillations computed in ref. (see figure 9a of that paper) with 5*17(0) = 
I9~^2^yh and the standard analysis procedure used here(cf figure ^ and in ref. p6[ . 
The right panel, which is the same as figure shows the allowed regions computed with 
the same procedure but with the Junghans et al. value of 517(0) = (22.3 it 0.9)eVb ||8|. 

The modification of the flux and uncertainties considered in this paper results 
in a reduction of the allowed neutrino oscillation parameter space. Since no pure 
sterile oscillation solutions are found at the 3a level, we focus on oscillations involving 
only active neutrinos. We find allowed regions for the LMA, LOW, and VAC solutions 
with our standard analysis strategy. The allowed regions for the preferred LMA and 
LOW solutions both become smaller with the larger adopted 6*17(0), and the fit to the 
LOW solution becomes relatively less good. The SMA and Just-so^ solutions, which 
are marginally allowed in the standard 'Before' panel of figure |^, are not present in 
the 'After' panel at the 3-a level. 

The differences between the Before and After panels are mainly due to two effects, 
which are described below. 

(a) The increase of the ^B normalization results in a smaller value for Xmim which 
is found for the LMA best-fit point. With the new value of 6*17(0), we find 
X^(LMA)^in = 34.5 while before we had x^(^^^)mm = 35.3. Although the 
theoretical errors are reduced in the new solution, the global minimum (which 
lies within the boundaries of the LMA solution) is "deeper." This effect can be 
confirmed by comparing the results of figure [l|a and figure |l|c, which were 
computed with and without the BPOO constraint on the ^B neutrino flux. 
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The global minimum points are much closer together in neutrino parameter 
space than they were when the older value of 5*17(0) was used (cf. results 
given in ref. |]16|)- With the new *B normalization the LMA best-fit point 
is not only the best-fit point in oscillation parameter space but it also corre- 
sponds very closely to the predicted neutrino flux given in table |1|. We 
find that 0(^3; LMA best - fit) = 0.9990(^6; BPOO + New^B), where the coef- 
ficient 0.999 applies for case (c) (free ^B analysis). For our standard analysis, 
case (a ), the coefficient is 0.96^. Presumably, this precise agreement between 
the best-fit LMA neutrino flux and the BPOO + New '^B flux is accidental, since 
the allowed Sa range for the ^B flux, within the domain of the LMA oscillation 
parameters, is (0.54 - 1.36)0(^B; BPOO + New^B). 

(b) The theoretical error for the *B flux is reduced according to eq. (|2.1| ) from an 
average value of 19% with the Adelberger et al. value of S'i7(0) to 14.5% 
with the Junghans et al. value of S'i7(0). The reduction of this theoretical 
error pushes the previously marginally allowed SMA and Just So^ solutions (for 
these somewhat unfamiliar solutions see the original discovery papers ||31| and 
ref. [|13|) over the edge (beyond the 3a allowed region). As discussed in section 
5.3 of ref. [0, the error on the predicted ^B neutrino flux played a crucial role 
in previously allowing, with the standard analysis prescription, the marginal 
SMA and Just So^ solutions. 

We have performed a series of numerical experiments to determine which effect, 
either the increase in the predicted '^B flux or the decrease in the uncertainty in 
the predicted flux, is the dominant effect in reducing the allowed parameter space. 
We find that the increase in the predicted flux contributes the most to worsening 
the flt for the LOW solution while the decrease in the uncertainty in the predicted 
flux is dominant reason why the SMA and Just So^ regions no longer appear at 3a. 
The increase in the predicted flux and the decrease in the theoretical uncertainty 
contribute comparably to reducing the size of the LMA allowed region. 

3.2 The allowed range for the '^B neutrino flux 

The analysis described above yields best-estimates and 3a allowed ranges for the 
total ^B neutrino flux. Let /b be the inferred ^B neutrino flux produced in the sun 
in units of the best-estimate predicted BPOO neutrino flux, 

(Py -DjBPOO+NewSB 



^Thcre is a simple reason why the coeiBcient, defined in eq. (3.2) as fe, deviates shghtly more 
from unity when (/)(^B) is constrained by the standard solar model than when it is not constrained. 
The reason is that the theoretical uncertainty is given as a fraction of the total *B neutrino flux. 
For higher fluxes, the theoretical error is larger and the value of is reduced somewhat more. As 



can be seen from eq. (3.2), fe depends inversely upon the theoretical flux. Therefore, higher values 



of the flux, i. e., lower values of fe, are slightly favored. 



- 12 - 



Table ^ gives the best-estimate and allowed range of the solar *B neutrino flux. 
The values in brackets in table ^ were obtained by our standard analysis, case (a), 
while the values without brackets were derived allowing the neutrino flux to be a 
free parameter, i. e., case (c). The bracketed range is always slightly smaller than the 
unbracketed range, which just reflects the fact that for strategy (a) the penalizes 
large differences from the standard solar model value of the flux while the ^B 
flux is unconstrained for strategy (c). However, the two strategies, (a) and (c), yield 
similar allowed ranges for /b^. 

In all cases, the inferred range of /b allowed by the global oscillation analysis 
using data from all of the experiments is smaller than the allowed range obtained 
by the SNO collaboration |^ from comparing the SNO and Super-Kamiokande total 
rate measurements. In our notation, the SNO collaboration found /b = 0.92 ± 0.50, 
where the result is expressed in terms of the BPOO flux calculated with the Junghans 
et al. value of 5*17(0) (see table |I| and ref. 0) and we quote the 3cr allowed range. Our 
result is /b = 0.88 ± 0.48 [or, for strategy (a), /b = 0.86 ± 0.38]. Of course, the SNO 
determination is more direct. The procedure described here assumes the validity of 
the two-neutrino oscillation analysis and uses solar model fluxes and uncertainties 
(see especially discussion in section |7?T|) . 

Even when the ^B neutrino flux is unconstrained in the oscillation analysis, 
almost all of the allowed range of /b is within the 3(T uncertainty of the BPOO solar 
model predictions. Only for the SMA solution does the allowed range of /b not 
overlap with the 3(T predicted range of the BPOO solar model. 

The procedure described above assumes the validity of the two-neutrino oscilla- 
tion analysis for pure active oscillations. As shown in ref. 1^^, the allowed range of 
/b can be modified by allowing for the possibility of an admixture of sterile neutri- 
nos in the oscillations within the adiabatic regime. In order to test this possibility 
we have repeated our global analysis (c) with an arbitrary active-sterile admixture 
within the region 10-3 < AmVeV^ < 10"^ , 0.05 < tan^^ < 5. We find that the 
maximum allowed ^B neutrino flux is increased to /b < 2.9 (at 3(j for 3 d.o.f.) and 
corresponds to oscillations into a 25% active and 75% sterile state. However, these 
mixed scenarios including sterile neutrinos give a worse fit to the data than pure 
active oscillations. A large component of sterile neutrinos does better if one includes 
in the analysis only the total event rates, as in ref. For completeness, we briefly 
describe how this is possible. In order to fit the appearance of active neutrinos 

^For strategy (a), we first obtain the allowed region of oscillation parameters using the BPOO 
values for the neutrino flux and uncertainties and then for each pair of (Am^ , tan^ &) within 
the allowed region we obtain the "optimum" value of /b that best fits the data. The range given in 
Table ^ is the range of the optimum /b values obtained with this procedure for all the oscillation 
parameters within the Scr allowed region in Fig. 1(a). For strategy (c), we defined the allowed regions 
allowing the *B neutrino flux to be a free parameter, unconstrained by the BPOO predictions. Thus 
the range given in Table ^ is the range of the optimum /b obtained for all the oscillation parameters 
within the 3(t allowed region in Fig. 1(c). 



- 13 - 



Solut ion /b ( Best ) 



/b(3(J range) 



LMA 1.00 0.54- 1.36 [0.57- 1.24] 

LOW 0.73 0.54- 0.93(0.55 - 0.91] 

VAC 0.59 0.48 - 0.68 [0.48 - 0.66] 

SMA 0.44 0.40 - 0.50 \ no solution] 



Table 4: Allowed values for the '^B neutrino flux. The table gives for each allowed 
oscillation solution the best-fit values for /b and the 3a range of allowed values. The 
quantity /b is the '^B neutrino flux expressed in units of the best-estimate BPOO value for 
the ^B neutrino fiux, see eq. ( |3.2D . The values of /b given in brackets were obtained by our 
standard analysis procedure, (a); the values presented without brackets were obtained by 
analysis procedure (c), in which the ^B flux is not constrained by solar model predictions. 



at Super-Kamiokande versus SNO, one needs to decrease Pee ^oi B neutrinos |p2 
(which is compensated by the increase of /b ), while in order to keep the agreement 
with gallium and chlorine rates the survival probability at lower energies must not be 
substantially modified. Barger et al. have pointed out in ref. p5[ that this tuning 
can be achieved within the adiabatic regime. However, as discussed in ref. |jl6|, one 
has less freedom for tuning in the global analysis. Within the adiabatic regime, the 
tuning of the probability can only be achieved by lowering both the value of tan^ 6 
and Am^, intruding into the region of a predicted large day-night variation in the 
Super-Kamiokande experiment. 

3.3 Methods of analysis 

In this subsection, we describe the implementations of the three different analysis 
strategies that were used to construct figure |l| and figure ||. For most readers, the 
outline of the analysis strategies given in section p.l.l| will be sufficient. Section ^ 
is intended only for aficionados of solar neutrino oscillation studies. 

The results shown in figure |l| and figure ^ were derived using the CC event rate 
measured by SNO 0] , the Chlorine and GaUium event rates |^ |, |^ (we use here 
the weighted averaged of the GALLEX/GNO and SAGE rates), and the 2 x 19 bins 
of the (1258 day) Super-Kamiokande |^ electron recoil energy spectrum measured 
separately during the day and night periods. We include, as described below, the 
predicted neutrino fluxes and their uncertainties, table ^ and table || of the present 
paper, for the standard solar model [|l^] (BPOO + New ^B). We use the distribution 



for neutrino production fractions and the solar matter density given in ref. |T2[ and 
tabulated in |http://www.sns.ias.edu/^jnb . 



In order to explore the robustness of the inferred allowed regions in neutrino 
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parameter space, we obtain the permitted regions using calculations based upon 
three different plausible prescriptions for the statistical analysis that have been used 
in the published literature and which we label as follows. 

(a) For our "standard" analysis, we adopt the prescription described in refs. [0, |T^]. 
We do not include here the Super-Kamiokande total rate, since to a large extent 
the total rate is represented by the flux in each of the spectral energy bins. We 
define the function for the global analysis as: 

Xglobal,a = X] {Rf - RT^)(^G%i^f ~ ^T"^)^ (3-3) 
ij=l,41 

where Gq^^ = aj^^j + crgp^y Here an^ij is the corresponding 41 x 41 error 
matrix containing the theoretical as well as the experimental statistical and 
systematic uncorrelated errors for the 41 rates while (Tsp,ij contains the assumed 
fully-correlated systematic errors for the 38 x 38 submatrix corresponding to 
the Super-Kamiokande day-night spectrum data. We include here the energy 
independent systematic error which is usually quoted as part of the systematic 
error of the total rate. The error matrix (TR^ij includes important correlations 
arising from the theoretical errors of the solar neutrino fluxes, or equivalently 
of the solar model parameters. We follow the principles described by the Bari 
group ||3^ in including the solar model uncertainties, using the updated data 
presented in table ^. 



(b) For the "all rates" analysis, we follow the prescription described in refs. [p2 



p3| , p4| , |T^, |3^. In this approach, we use all the total event rates in the anal- 
ysis, including the Super-Kamiokande total event rate with its corresponding 
theoretical error. We define the function for this global analysis as: 



Xglobal,b — XRates.b + XsP,b5 (3-4) 



where 



XRates,b 



J2 (^f - l^n^'n^^.iRf - ^r)' (3-5) 

xip.b = E («spi?f - i?r)^5P..(«sp^f - Rfn- (3.6) 

i,j=l,38 

In XRates5 include the CC event rate measured at SNO 0, the Chlorine 
and Gallium p, ^ Q event rates and the Super-Kamiokande total event rate. 
The matrix aji^ij contains the theoretical uncertainties, as well as the experi- 
mental statistical and systematic errors, for the total rates. In particular, (TH,ij 
includes the energy independent systematic error for the Super-Kamiokande 
rate. On the other hand, the matrix crsp,jj contains both the statistical as well 
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as the systematic errors (both those that are correlated with energy and those 
that are uncorrelated) of the 38 spectral energy data; asp^ij does not include 
the theoretical error of the flux. 

In order to avoid double counting, we allow a free normalization (for which no 
theoretical error is included) for the Super-Kamiokande day and night electron 
recoil energy spectra. Following the usual procedure for this approach, we in- 
troduce an arbitrary normalization factor, asp, for the amplitude of the energy 
spectra. For each value of the oscillation parameters Xgiobai,6 minimized with 
respect to asp- 

(c) The "free ^B" analysis is identical to our standard analysis, described in (a) 
above, except that the ®B neutrino flux is not constrained by solar model pre- 
dictions. This freedom is implemented by multiplying the ^B flux contribution 
to the Rj^ in eq. ( |3.3| ) by the normalization factor, /b, and removing the theo- 
retical ^B flux errors from cr/j,jj. In this case, the resulting function can also 
be written as 

Xglobal.c — Xptates.c + XsP.C (3-7) 

where 

xLtcc = E «(/b) - ^r)^i^i(^f (/b) - i?r)' (3-8) 

Xip,c = E (^*'(/b) - i?r)^5P..(^f (/b) - ^D- (3-9) 

jj=l,38 

In Xptates c include only the CC event rates measured in the SNO, Chlorine, 
and Gallium experiments. For these three rates, the error matrix crn^ij only 
differs from ct/j y due to the absence of the theoretical error for the ^B neutrino 
flux. The error matrix d'sp,ij only differs from o'sp^ij by the inclusion of the 
energy independent systematic error, which is usually quoted as part of the 
systematic error of the total rate. 

For each value of the oscillation parameters, Xgiobai b is minimized with respect 
to /b. Notice that the only difference in the "free ^B" analysis when performed 
in terms of BPOO or in terms of BPOO + New ^B is a shift in the corresponding 
best value of /b for each value of the oscillation parameters /b — > /b (fff); the 
allowed regions are left unchanged as can be seen by comparing Fig. |l]c with 



the right panel of Fig.9 in Ref. |]T6|. 



4. Predictions for SNO 

In this section, we present the predictions for the Sudbury Neutrino Experiment 
(SNO), ref. [0, of the currently allowed neutrino oscillation solutions (cf. flgure |^) 
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for the neutral current to charged current ratio (section [4.1|) , for the CC day-night 
effect (section [4.2|) , for the z/g + e day-night effect (section ^l3|) , and for the first and 
second moments of the electron recoil energy spectrum (section |4.4[ ). In Section 



we show that the neutrino oscillation parameters will be well determined if we are 
lucky and SNO measures a large value (> 3.3) of the neutral current to charged 
current ratio or a large value (> 0.1) for the day-night effect. 

The SNO collaboration is studying charged current (CC) neutrino absorption by 
deuterium, 

Ve + d^p + p + e~ , (4.1) 
as well as the neutral current (NC) neutrino disassociation of deuterium, 

Vx + d n + p + ■u'.j. , (x = e,yU, r). (4.2) 

SNO is the only operating (or completed) solar neutrino experiment that determines 
a CC rate for electrons whose energies are measured. The SNO neutral current 
detection is also unique. In addition, SNO measures neutrino-electron scattering, 

z/a; e ^ e + i/^. , {x = e,n, t). (4.3) 

The statistical uncertainties for neutrino-electron scattering in the SNO detector are 
less favorable than those available from the Super-Kamiokande experiment [Q. 

The predicted measurable effects for SNO have been investigated in detail in 
ref. for the larger range of neutrino oscillation solutions that were viable prior to 
the publication of the first SNO results on the rate of the CC reaction. We follow 
here the notation and the techniques described in ref. |3^, concentrating on those 



effects that were shown previously to be the most likely to be measurable. We do 



not repeat here calculations of the effects that were discussed in ref. |^ and found 
to be small, such as the neutral current day-night effect or the seasonal difference 
(winter-summer difference). 

For the CC reaction, we adopt a 6.75 MeV kinetic energy threshold for the recoil 
electrons, the threshold used by the SNO collaboration in their CC measurement 
reported in ref. 0, as the standard value used in the figures. In order to illustrate 
the dependence upon the CC threshold, we also present results in the tables that 
refer to a 4.5 MeV recoil electron kinetic energy threshold. 

4.1 The Neutral Current to Charged Current Double Ratio 

In this subsection, we present predictions for the ratio of neutral current events (NC) 
to charged current events (CC) in SNO. The most convenient form in which to discuss 
this quantity is PB| 



[NC] _ ((NC)obs/(NC) 
[CC] ((CC)obs/(CC) 
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The ratio [NC]/[CC] is equal to unity if nothing happens to the neutrinos after 
they are produced in the center of the sun (no oscillations occur). Also, [NC]/[CC] 
is independent of all solar model considerations provided that only one neutrino 
source, ^B, contributes significantly to the measured rates. Inserting into eq. (25) 
of ref. |]35| the 3a upper limit measured for the hep flux by the Super-Kamiokande 
collaboration |^, one finds that hep neutrinos affect the value of [NC]/[CC] by less 
than 0.05% for a 4.5 MeV CC kinetic energy threshold (and by less than 2% for 
a 6.75 MeV CC threshold.) Moreover, the calculational uncertainties due to the 
interaction cross sections and to the shape of the neutrino energy spectrum are 
greatly reduced by forming the double ratio |^ ^ . 

Figure ^ and table || present the calculated range of the double ratios for the 
oscillation solutions that are currently allowed at 3a. The table gives the best-fit 
values for [NC]/[CC] as well as the maximum and minimum allowed double ratios 
for a threshold recoil electron kinetic energy for the CC reaction of 4.5 MeV and 
separately for a CC threshold of 6.75 MeV. The calculated double ratio is insensitive 
to the CC threshold within the range of thresholds that are considered. 

The la predicted range for our standard global analysis strategy (a) is 3.45lo;54 
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Figure 3: The neutral current to charged current double ratio, [NC]/[CC] . The 
double ratio, [NC]/[CC], is defined by eq. (4^). The three panels of predictions in figure ^ 
were derived from the 3a ranges shown in the three panels of global solutions illustrated in 
figure H and obtained using the three analysis strategies described in section Figure |^ was 
constructed using a 6.75 MeV total electron energy for the CC threshold. Table ^ shows 
that, within the allowed domain of neutrino oscillation solutions, the value of [NC]/[CC] is 
not sensitive to the threshold energy adopted for the CC reactions. The tiny 3a estimated 
measuring error in the left hand panel assumes the systematic uncertainties projected in 
ref. 1 35 1 and 3000 CC events and 1200 NC events. The ultimate SNO measurement will 
probably be even more precise than this estimate. 
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b.f. max mill b.f. max min 

Scenario 4.5 MeV 4.5 MeV 4.5 MeV 6.75 MeV 6.75 MeV 6.75 MeV 





[NC]/[CC] 


[NC]/[CC] 


[NC]/[CC] 


[NC]/[CC] 


[NC]/[CC] 


[NC]/[CC] 


LMA 


3.44 


5.28 


1.79 


3.45 


5.28 


1.82 


LOW 


2.39 


3.22 


1.71 


2.37 


3.20 


1.71 


VAC 


1.76 


2.06 


1.43 


1.82 


2.17 


1.46 




An-d 


An-d 


An-d 


An-d 


An-d 


An-d 


LMA 


7.4 


19.5 


0.0 


8.3 


21.4 


0.0 


LOW 


4.3 


10.4 


0.0 


3.7 


9.5 


0.0 


VAC 


0.1 


0.3 


-0.2 


0.2 


0.5 


-0.3 



Table 5: SNO Neutral Current to Charged Current Double Ratio and Day- 
Night CC Asymmetry. The table presents the double ratio, [NC]/[CC] and ^n-d 
(in %). The results are tabulated for different neutrino oscillation scenarios and for two 
different thresholds of the recoil electron kinetic energy used in computing the CC ratio, 
4.5 MeV (columns two through four) and 6.75 MeV (columns five through seven). The 
ranges are obtained for the 3cr regions for analysis (a). 



for a CC threshold of 6.75 MeV(3.44l°:^° for a 4.75 MeV threshold). The predicted 
range for analysis strategy (b) is slightly smaller and for analysis strategy (c) is 
somewhat larger ( 3.65lQ;go). 

The best-fit values for the double ratio for oscillations into active neutrinos range 
between 1.8 and 3.4 (1.8 and 3.5) for a 4.5 MeV (6.75 MeV) CC threshold. The 
maximum predicted value for [NC]/[CC] is an enormous 5.3 for an extreme LMA 
solution. The minimum calculated value for the double ratio is 1.4, which is implied 
by an extreme VAC solution. 

Figure ^ compares the "Before" and "After" predicted values for the neutral 
current to charged current double ratio. For the favored three solutions, LMA, 
LOW, and VAC, the predicted double ratios are not very sensitive to the precise 
value of 517(0) that is adopted. We note however that in all three cases adopting the 
Junghans et al. value of 517(0) causes the values of [NC]/[CC] to become slightly 
larger (more distant from the no oscillation value of 1.0) as a consequence of the 
reductions and shifts of the allowed regions. 

4.2 The Charged Current Day-Night Effect 

In this subsection, we discuss the difference between the charged current event rate 
observed at night and the charged current event rated observed during the day. This 
difference in event rates has been evaluated previously for a variety of experiments 
by many authors, including those listed in refs BO, O, pi], . 
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Figure 4: [NC]/[CC] "Before and After." The left panel shows the 3a" allowed ranges 
for the neutral current to charged current double ratio using the neutrino oscillations 
solutions computed in ref ||l^ with 517(0) = lOlgeVb |19| and the standard analysis 



procedure used here(cf. figure IT]) and in ref. [16|. The right panel, which is the same as 



figure shows the allowed values for [NC]/[CC] computed with the same procedure but 
with the Junghans et al. value of ^^(O) = (22.3 ± 0.9)eVb @. 

We concentrate here on the difference, An-d, between the nighttime and the 
daytime CC rates for SNO, averaged over one year. The definition of v4n-d is 

[Night + Day] ^ ' 

The value of An-d is zero for neutral current detection of oscillations into active 
neutrinos. (For neutrino oscillations into sterile neutrinos, An-d has been calculated 
in ref. and shown to be small, less than 1% for the range of neutrino parameters 
allowed before the SNO experimental results were available.) 

Table |^ and figure |^ present the range of predicted differences between the aver- 
age rate at night and the average rate during the day [i.e., 100 x y4N-D(SN0 CC) of 
eq. (|4.5| )]. The calculated predictions in table ^ are given for a 4.5 MeV and an 6.75 



MeV CC electron recoil kinetic energy threshold. 

For most of the MSW oscillation solutions, the predicted day-night differences 
are only of order a few percent. However, for the LMA solution, the predicted 
difference can reach as high as 20% for a 4.5 MeV threshold (21% for a 6.75 MeV 
threshold). For very special choices of the LOW parameters, An-dISNO CC) can be 
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Figure 5: The percentage difference between the night and the day CC rates. 

The figure shows the night-day percentage difference for the charged current rate in SNO, 
i.e., 100 X ^N-D defined in eq. (4.5). Predictions are shown for the solar neutrino oscillation 



scenarios allowed at 3a and illustrated in figure ^. The three panels refer to results for 
different analysis strategies described in section ^. Figure |5| was constructed using a recoil 
electron kinetic energy threshold of 6.75 MeV for the CC events. Table |5| gives numerical 
results for ^n-d for kinetic energy thresholds of 4.5 MeV and 6.75 MeV. 

as large as 9% or 10%. For vacuum oscillations, there is a small day-night effect that 
is due to the dependence of the survival probability upon the earth-sun distance. 
This dependence has been calculated in ref. |^5| and corresponds in all the allowed 
cases to |An-d| < 1%- 

The la predicted range for our standard global analysis strategy (a) is An-d 
(SNO CC) = 8.3^^:° % for a 6.75 MeV CC threshold and An-d(SNO CC) = 7Attl % 
for a 4.5 MeV threshold. The predictions for analysis strategies (b) and (c) are very 
similar. 

Figure P shows the "Before and After" comparison between the range of predicted 
values for 74n-d(SNO CC) using the older S'i7(0) and the most recent and accurate 
determination of 6*17(0). For the three favored solutions, LMA, LOW, and VAC, 
the allowed range of ylN--D(SNO CC) is only changed slightly, the maximum value is 
increased (decreased) a small amount for the LMA (LOW) solution. 

4.3 The u + e Scattering Day-Night Effect 



We present in section |4.3.1| the calculated day-night effects for neutrino-electron 



scattering in the SNO detector that are predicted by the current globally-allowed 



oscillation solutions. In section |4.3.2| , we present and discuss the predicted corre- 
lation, which depends somewhat on the particular oscillation solution, between the 
day-night effect in the CC event rate and the day-night effect in the neutrino-electron 
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Figure 6: Night-Day CC "Before and After." The left panel shows the 3a allowed 
ranges in percent for the average CC night minus the day event rates, ^n-d, as defined in 
eq. ( [4.5p ; the values in the left panel were computed using the neutrino oscillations solutions 
in ref |]l6[ obtained with S'i7(0) = 19l26Vb and the standard analysis procedure used 
here(cf. figure 0) and in ref. The right panel, which is the same as figure |5^, shows 
the allowed values for j4n-d computed with the same procedure but with the Junghans et 
al. value of S'i7(0) = (22.3 ± 0.9)eVb g. 



scattering event rate. This correlation constitutes an important consistency check 



on the oscillation solution [43 



4.3.1 Predicted values for v4n_d(SN0 ES) 

Figure |^ and table H show the predicted (percentage) night-day differences foiu + e 
scattering in the SNO detector. The calculated values for An-d(SNO ES) were 
obtained using the 3a boundaries of the global oscillation solutions shown in figure ^ 
The values of An-d( ES) that are predicted for the Super-Kamiokande and the 
SNO locations are very similar. For the Super-Kamiokande location, the LMA 
solution, and a kinetic energy threshold of 4.5 MeV ( corresponding to the re- 
cent Super-Kamiokande electron recoil energy threshold 0), the predicted range 
of An-d(SK ES) is between 0.0% and 10.4%, with a best-fit value of 4.3% . Simi- 
larly for the LOW solution, the predicted range is between 0.0% and 7.3%, with a 
best-fit value of 3.1%. The predicted values of y4N_D( ES) at both the SNO and the 
Super-Kamiokande locations are negligibly small for the VAC solution. Comparing 
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Figure 7: The percentage difference between the night and the day z/+e scatter- 
ing rates. The figure shows the night-day percentage difference, i.e., 100 x ^n-d(SNO ES) 
defined in eq. ( |4.5| ). Predictions are shown for the solar neutrino oscillation scenarios al- 
lowed at 3(7 and illustrated in figure |^. The three panels refer to results for different 
analysis strategies described in section ^. Figure |5| was constructed using a recoil electron 
kinetic energy threshold of 6.75 MeV. Table ^ gives numerical results for An-d(SNO ES) 
for kinetic energy thresholds of 4.5 MeV and 6.75 MeV. 

with table we see that the range of predicted values of v4n-d( ES) differs, between 
the SNO and Super-Kamiokande locations, by at most 7% of the predicted range. 

The predictions for SNO can be compared with the measured result reported by 
the Super-Kamiokande collaboration 

^n-d(SKES) = 0.033 ±0.022(stat.)+°:°}^(sys.). (4.6) 
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max 


min 


b.f. 


max 


min 


Scenario 


4.5 MeV 


4.5 MeV 


4.5 MeV 


6.75 MeV 


6.75 MeV 


6.75 MeV 




An-d 


An-d 


An-d 


An-d 


An-d 


An-d 


LMA 


4.1 


10.1 


0.0 


4.7 


11.4 


0.0 


LOW 


3.3 


7.8 


0.0 


2.9 


7.1 


0.0 


VAC 


0.0 


0.1 


-0.1 


0.1 


0.3 


-0.2 



Table 6: The predicted u + e Day-Night Asymmetry in SNO. The table presents 
^n-d(SNO ES) in percent. The results are tabulated for different neutrino oscillation 
scenarios determined by our standard strategy and for two different thresholds of the recoil 
electron kinetic energy, 4.5 MeV (columns two through four) and 6.75 MeV (columns five 
through seven) . The ranges are obtained for the 3a regions for figure ||(a) . 
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Not surprisingly, the best-fit predictions for tlie u+e day-niglit scattering difference in 
tlie SNO detector implied by the global solutions shown in figure |l|a are all rather close 
to the best-estimate value obtained by Super-Kamiokande. The SNO predictions are 
also bounded by the 3a Super-Kamiokande limits, i.e., the predicted SNO values 
lie between 0% and 11%. Figure |^ shows that the LMA solution allows the largest 
values for An-d(SNO ES) and the VAC solution predicts no day-night difference in 
u + e event rates. 

The la allowed range for v4n-d(SN0 ES) is 4.7l|I % for 6.75 MeV and 4.1+|^ % 
for a 4.5 MeV threshold, both for analysis strategy (a). The results for strategies (b) 
and (c) are very similar. 

The "Before and After" predictions for 74n-d(SNO ES) are essentially the same, 
with only minor differences. That is to say, that the results do not depend signif- 
icantly on the choice between the previously standard value of 5*17(0) [|1^] and the 
most recent and precise value of 5*17(0) 0. 

4.3.2 Correlation between An-d(SNO CC) and An d(SNO ES) 

Figure ^ shows the correlation between the predicted values for the day-night effect 
of the CC rate, An-d(SNO CC), and the day-night effect of the neutrino-electron 
scattering rate, 74n-d(SNO ES). The two day-night effects are essentially propor- 
tional to each other. The essential features of the correlation shown in figure |^ can 
be derived analytically, as we show in the following discussion. 

The survival probability of electron neutrinos can be written as 

Pee = Pd - (1 - 2P,) • COS 203 " /reg, (4.7) 

where 

Pb = - + -(l-2Pc) ■cos2^5-cos2^ (4.8) 

is the survival probability in the absence of the Earth-matter effect, i.e., during 
the day. Here 6s is the matter mixing angle at the production point inside the Sun 
{9s = f when the production point occurs at much higher densities than the resonant 
point) , and Pc is the level crossing probability which describes the non-adiabaticity 
of the conversion inside the Sun. The Earth regeneration factor, /reg, is defined as: 

/reg = P2e ' sill' 9 , (4.9) 

where P2e is the probability of the 1^2 — * i^e conversion inside the Earth. In the 
absence of the Earth-matter effect we have P2e = sin^ 9, so that /reg = 0. 

For both LMA and most of LOW region the adiabaticity condition is satisfied 
and we can take P^ ~ 0. Furthermore in the LOW region and in the lower part of 
the LMA region one has cos 2^5 = — 1. With these approximations the probabilities 
take the form 

Pee = Sin2^ + /reg , (4.10) 

PD = sm'^9, (4.11) 
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Figure 8: The correlation between the day-night effects for v + e scattering and 
for CC events. The day-night effects shown in figure |5| and figure ^ are correlated for the 
standard analysis strategy that resulted in figure ||a. The predictions for the VAC solutions 
are essentially a point at the position (0,0). 



And the day-night asymmetry for CC events can be written as 

(/rcg) 



An-d(SNO CC) = 2^^^— ^ = 2 

Vee) 



sin e + (/reg) 



(4.12) 



where we have denoted by () the averaging over the neutrino spectrum, the interac- 
tion cross section, and the energy resolution. 

The [ES] event rate is [ES] ~ (Pee + '"(I — Pee)) where r = a^/ae — 0.15 is the 
ratio of the the t'e — e and v^ — e elastic scattering cross-sections. Hence the day- night 
asymmetry for neutrino-electron scattering events can be written as 



An-d(SNO ES) = 2 



^)(/rog) 



sin 9 + r cos^ 9 



(4.13) 



From eq. ( [4.12[ ) and eq. ( [4.13| ), we see that the slope in the y4N-D(SN0 CC) versus 
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An-d(SNO ES) plot is given by 

ylN-D(SNO ES) ^ sin^g+(/,eg) sin^g 

An-d(SNOCC) ^ ^^sin2^ + rcos2e + (1 -r)(/reg) ^ ^ sin^ ^ + r cos2 ^ 

(4.14) 

where the last implication applies in the region of small asymmetries. Using this 
relation and the ranges of mixing angles in the LOW and LMA solution we can 
reproduce the slope of the regions. From eq. ( |4.14| ), we see that the ratio of the 



asymmetries increases with the mixing angle and is larger for mixing angles on the 



dark side {6 > 7r/4). Therefore eq. ( 4.14 ) explains why the LOW solution has a 



larger slope in figure |^ than the LMA solution. 

4.4 CC recoil energy spectrum: first and second moments 

The solar influence on the shape of the '^B neutrino energy spectrum is only about 
one part in 10^ MM. Therefore, in the absence of neutrino oscillations or other 



new physics, the energy spectrum of solar neutrinos incident on terrestrial detectors 
should be the same as the energy spectrum of *B neutrinos produced in the labora- 
tory. The incident solar neutrino energy spectrum can be studied well by observing 
the recoil energy spectrum of electrons produced by neutrino absorption (CC) re- 
actions on deuterium, as shown in eq. (|4.1|) . In the CC reaction, nearly all of the 
available kinetic energy is taken by the recoiling electron since the other particles in 
the final state are (massive) baryons. 
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Figure 9: The first moment: the fractional shift, ST, in percent for the average 
electron recoil energy. The figure shows the fractional change in percent for the average 
electron recoil energy, (T), relative to what is produced by an undistorted neutrino 
energy spectrum. The predictions are shown at 3a C.L. for the three different analysis 
strategies for neutrino oscillation solutions that are illustrated in figure |^. The calculations 
were performed assuming a 6.75 MeV recoil electron kinetic energy threshold. 
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Figure 10: The second moment: the fractional shift, (5cr, in percent for the 
standard deviation of the electron recoil energy spectrum. The figure shows the 
fractional change in the standard deviation of the electron recoil energy spectrum relative 
to what is expected for an undistorted neutrino energy spectrum. The predictions are 
shown at 3a C.L. for the three different analysis strategies for neutrino oscillation solutions 
that are illustrated in figure ||. The calculations were performed assuming a 6.75 MeV recoil 
electron kinetic energy threshold. 



The Super-Kamiokande experiment has already provided beautiful data on the 
recoil energy spectrum from neutrino electron scattering . In recent publications by 
the Super-Kamiokande collaboration, these data are divided into 19 separate energy 
bins. A number of authors have questioned whether using so many energy bins in 
a analysis of all the available data gives too much weight to the measurement of 
the spectral parameters (see, e. g., ref. for a particularly insightful discussion of 
the statistical questions). 

In this subsection, we consider an alternative method of analysis of the SNO 
recoil energy spectrum in terms of small number of the lowest order moments of the 
distribution, a method that has been discussed in connection with solar neutrino 



experiments in refs. p5| , p^ , We follow the notation and analysis of refs. 



Figure ^ and figure |T0| show the computed first and second moments of the recoil 
electron energy spectrum from the CC reaction, eq. ( |4.1|) , for the three different 
analysis strategies discussed in section |^ and illustrated in figure |1[ These moments 
may well represent the essential physical content of the electron energy spectrum 
(cf. refs. ^) and could represent the results of the spectral measurements in a 



global oscillation analysis. 

For both the first and second moments, the VAC solution predicts the largest 
deviation from the undistorted spectrum. Nevertheless, the range of predicted values 
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is only about ±2% for the first moment and ±7% for the second moment. The non- 
statistical uncertainties in measuring the first and second moments in SNO have been 
estimated, prior to the operation of the experiment, in ref. and are, respectively, 
about 1% and 2%. 

Assuming the approximate validity of the prior estimate of the measuring un- 
certainties given in ref. [^], it seems very unlikely that SNO will determine a 3a 
deviation from the undistorted spectrum if either the favored LMA or LOW oscilla- 
tion solutions shown in figure |l| is correct. Although this is a negative prediction, it 
is nevertheless an important prediction. Even for the VAC solution, it will be very 
difficult to measure a significant deviation from an undistorted spectrum. The best 
opportunity to see a deviation will be in the second moment. 



4.5 SNO: What if we are lucky? 

Figure ^, figure ^, and table ^ show that if Nature has arranged things favorably, 
then only a rather small range of neutrino oscillation parameters can correctly pre- 
dict the results. If SNO measures either a large value for the neutral current to 
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Figure 11: If nature is kind. The left panel shows the currently allowed regions for the 
neutrino oscillation parameters that predict a neutral current to charged current double 
ratio, [NC]/[CC], that is > 3.3 and the right panel shows the currently allowed region 
that predicts a CC night-day difference, An-d(SNO CC), that is > 0.1. The calculations 
were performed assuming a 6.75 MeV recoil electron kinetic energy threshold. The regions 
shown in figure 11 will change after measurements become available on either (or both) 
[NC]/[CC] or An-d(SNO CC) and a new global solution is calculated. 
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Figure 12: The correlation between the neutral current to charged current dou- 
ble ratio, [NC]/[CC], and the charged current day-night effect An dISNO CC). 
The figure shows the allowed values at 3cj of [NC]/[CC] versus 100 x An-d(SNO CC), 
using the analysis strategy (a) and a CC threshold of 6.75 MeV. The dotted lines show the 



limiting values ^n-d(SNO CC) = 10% and [NC]/[CC] = 3.3 used in figure 11 



charged current ratio, [NC]/[CC] > 3.3, or a large value for the night-day difference, 
An-d(SNO CC) > 0.1, then the Am^ and tan^^ will be rather well determined. 

Figure ^ shows the small range of neutrino parameters that, given the currently 
allowed regions defined by the currently available experimental data, can correctly 
predict the results if SNO shows that one of these inequalities is correct. In that 
fortunate situation, SNO will narrow down the allowed neutrino parameters to a small 
region in neutrino parameter space. With the currently available data, the boundaries 
of the region corresponding to [NC]/[CC] > 3.3 are for our standard analysis strategy, 
(a): 0.2 < tan^ 9 < 0.4 and 1.8 x lO'^ eV^ < Am^ < 1.0 x 10"^ eV^. For An_d > 0.1, 
the corresponding boundaries are 0.2 < tan^6' < 0.7 and 1.8 x 10"^ eV^ < Am^ < 
3.5 X 10^^ eV^. The only change in the above numbers if we use the free '^B strategy 
(cf. figure |T]c), rather than our standard analysis strategy is that the v4n_d range is 
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slightly affected and the upper limit for Am^ decreases to 3.3 x 10 ^ eV^. 

The regions corresponding to [NC]/[CC] > 3.3 and An-d(SNO CC) > 0.1 de- 
pend upon the data available when a global analysis is made. The regions shown 
in figure |TT] will change after measurements become available on either (or both) 
[NC]/[CC] or An-d(SNO CC) and a new global solution is carried out. 

If bi-maximal mixing |3^ is correct, one would not expect either [NC]/[CC] > 3.3 
(this would be a lOcr deviation from predictions based upon existing experimental 
data, cf. figure |l|) or y4N-D(SN0 CC) > 0.1 (a 3.8cr [3.1cr]deviation from the predic- 
tions of currently allowed solutions for analysis (a) [(c)]). 

Figure |12| shows the correlation between the currently allowed values of [NC]/[CC] 
and ^n-d(SNO CC). There is a general tendency that larger values of v4n_d 
(SNO CC) are associated with larger values of [NC]/[CC], although the correlation is 
rather broad. For a given value of [NC]/[CC], a large range of values ofy4N-D(SNO CC) 
is currently allowed and vice- versa. Figure ^ also shows clearly that for smaller val- 
ues of [NC]/[CC] and An-d(SNO CC) discriminating among different oscillation sce- 
narios will not be easy. For example, for the LMA solution, [NC]/[CC] < 2.1 implies 
y4N-D(SN0 CC) < 1% . On the other hand, there is a range of neutrino parameters 
characterizing the LOW solution for which one expects [NC]/[CC] < 2.1 while there 
exists a small but non vanishing day-night asymmetry 1% < y4N-D(SN0 CC) < 8%. 



5. Predictions for ^Be rate and day-night effect 

In this section, we present the predictions of the globally allowed solutions for the 



neutrino electron scattering rate and the day-night effect in the BOREXINO |]TU| 



and KamLAND [|TT| detectors. For specificity, we present results for the location 
of BOREXINO, but comment in the text on the differences in the predictions for 
KamLAND and BOREXINO. We assume the same energy window for recoil electrons 
for both detectors, so the only cause for the differences in the predictions is the 
different latitude of the detector locations. In all cases, the differences are smaller 
than the expected measurement uncertainties. 

We consider only electrons with recoil kinetic energies in the experimentally 
preferred range of 0.25 MeV < Tc < 0.8 MeV (cf. ref. [0]). For simplicity in presen- 
tation and following the general practice in the literature, we refer to these events 
as ''Be neutrino events. However, we include in our discussion neutrinos from all 
solar neutrino sources that produce recoil electrons in the indicated energy range. 
We make use of the predicted neutrino fluxes and uncertainties given by the BPOO 
standard solar model [|l^]. (The ^B neutrinos do not contribute significantly in this 
energy range.) The ^^N and ^^O neutrinos are, after ^Be, the next most impor- 
tant neutrino sources for producing electrons with recoil kinetic energies between 
0.25 MeV and 0.8 MeV. 
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b.i. 


max 


min 


b.i. 


max 


min 


Scenario 


[R(^Be)] 


[R(^Be)] 


[R(^Be)] 


An-d 


An-d 


An-d 


LMA 


0.65(0.65) 


0.76(0.76) 


0.58(0.58) 


0.0(0.0) 


0.1(0.1) 


0.0(0.0) 


LOW 


0.62(0.62) 


0.74(0.74) 


0.54(0.54) 


27(27) 


42(42) 


0.0(0.0) 


VAC 


0.63(0.62) 


0.79(0.78) 


0.53(0.51) 


-3.5(-3.8) 


-1.1(-1.3) 


-4.8(-5.3) 



Table 7: ^Be: Neutrino-electron scattering rate and Day-Night Asymmetry. 

The table presents the reduced neutrino-electron scattering rate, [^Be] (eq. p| ), and the 
night-day difference, An-d (see, eq. 4.5 but the result is presented here in %), for recoil 
electrons with kinetic energies in the range 0.25 MeV < Tc < 0.8 MeV. The numbers in 
parentheses were calculated assuming that only ^Be neutrinos contribute while the entries 
not in parentheses were obtained assuming the correctness of the BPOO standard solar 
model neutrino fluxes. The results are given for the three neutrino oscillation scenarios 
described in figure |a. The ranges shown in the table correspond to the 3a allowed regions 
for figure |l]a. 



In order to indicate the insensitivity of our conclusions to the poorly known 
CNO neutrino fluxes, table ^ gives results for the scattering rate and the day-night 
difference both with and without including CNO and pep neutrinos. The predicted 
values are insensitive to the non-'^Be neutrinos, if their fluxes are comparable to the 
BPOO standard solar model fluxes. 

We begin by discussing the u — e scattering rate in section lO and then discuss 



the day-night effect in section |5.2| . The hypothetical accuracy with which KamLAND 
and BOREXINO may, together with previous solar neutrino experiments, determine 
neutrino oscillation parameters has been discussed in ref. EH . 



5.1 ^Be V — e scattering rate 



Figure |I3| shows, for the currently allowed oscillation solutions, the expected 3a range 
in BOREXINO of the reduced neutrino-electron scattering rate, [''Be], 



'Be] 



(Standard Model Value)i 



0.25 MeV < Tc < 0.8 MeV. 



(5.1) 

The neutrino-electron scattering rate, fBe], is deflned in eq. ( ^.1|) as the ratio of 
the number of recoil electrons with kinetic energies in the range 0.25 MeV < Tg < 
0.8 MeV divided by the number that is expected if the BPOO model is correct and 
there are no neutrino oscillations. The summation over the different neutrino sources 
i includes all solar neutrino fluxes, although the [^Be] neutrinos are expected to 



- 31 - 



0,9 
0.8 
0,7 



^0-6 
0) 
CD 



0.4 
0.3 
0.2 



(a) 



(b) 



ic) 



LMA LOW VAC 



LMA LOW VAC 



LMA LOW VAC 



SMA 



Neutrino Scenario 



Figure 13: The '^Be neutrino event rate. The figure shows ['''Be], the calculated 
I'e + e~ event rate for all recoil electrons (including CNO neutrinos, see text) with kinetic 
energies in the range 0.25 MeV < Te < 0.8 MeV relative to the rate that is expected for 
the BPOO standard model neutrino flux and no neutrino oscillations. The three panels of 



predictions in figure 13 were derived from the three panels of global solutions illustrated in 
figure H using the three analysis strategies described in section ^. 



dominate according to the BPOO solar model. The n and r neutrinos are denoted by 
in eq. (§]T|). 

Figure |13| shows that all three of the analysis strategies described in section || 
yield rather similar predictions for the scattering rate. 

The predicted 3a range of the neutrino electron scattering rate is (cf. table 0) 
[^Be] = 0.65ltai27 where both the minimum value and the maximum value are 
achieved by vacuum solutions and the quoted limits correspond to our standard 
analysis strategy (cf. figure |l|a). For the favored LMA solution, the range is some- 
what smaller: [^Be] = 0.65tom (see also ref. [|^. The SMA solution, which is allowed 
at 3cr only in the "free ^B analysis strategy (cf. figure |l|c), predicts much smaller 
values for fBe], namely, [^Be] < 0.34 . 

The la predicted range is [^Be] = 0.65^q;q3 for our standard global analysis 
strategy (a). For the free *B analysis strategy (c), the spread in predicted values 
of [''Be] is somewhat larger, [''Be] = 0.66]'^o;q4. Analysis strategy (a) and (b) predict 
essentially the same la range. 

For KamLAND, the predicted range of [^Be] is essentially the same as for 
BOREXINO. For the LMA and VAC solutions, the BOREXINO and KamLAND 
best-fit and maximum and minimum values are the same to three significant figures. 
For the LOW solution, the predicted range for KamLAND is about 1% lower than 
for BOREXINO; the best-fit, maximum and minimum values are shifted downward 
for KamLAND. 
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Figure 14: The percentage difference between the night and the day rates for 
recoil electrons with kinetic energies in the range 0.25 MeV < Tg < 0.8 MeV. 

The figure shows the night-day percentage difference, i.e., 100 x ^n-d defined in eq. (^). 
The predictions are given for the solar neutrino oscillation scenarios allowed at 3a and 
illustrated in figure |^. The three panels refer to results for different analysis strategies 
described in section ^. 



5.2 ^Be day-night variations 



Figure 14 shows the predicted percentage difference between the night and the day 



rates that is expected for the BOREXINO experiment. As a number of previous 



authors have discussed (see, e. g., refs. [|T^, [5^, ^ and references cited therein), the 
LOW solution is the only currently allowed oscillation solution that predicts a large 
night-day difference in the rates. The LMA solution predicts a negligible variation 
and therefore the la predicted range is effectively zero. 

The small night-day difference for the VAC solution shown in figure n and 
table 1^ is the first calculation of this variation, with which we are familiar, for the 
low energies to which BOREXINO and KamLAND are sensitive. The physical origin 



of this variation has been described in section IXA of ref . fS^j ; it is due to the fact that 
the VAC survival probabihty depends upon the earth-sun separation and the longest 
nights occur in the Northern hemisphere when the earth-sun distance is shortest. 

For KamLAND, the predicted range of values of v4n-d is very similar to the range 
predicted for BOREXINO. For the LMA solution, the predicted range of An_d for 
the two experiments is the same to an accuracy of better than 0.1%. For the LOW 
solution, the maximum value for KamLAND is about 2% less than for BOREXINO 
(41% compared to 42%) and for the VAC solution, the minimum value is —3.9% for 
KamLAND compared to -4.8% for BOREXINO. 
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6. Predictions for KamLAND reactor experiment 



The KamLAND detector |TT|, located in the site of the famous Kamiokande 
experiment |Q , consists of approximately one kton of liquid scintillator surrounded 



by photomultiplier tubes. KamLAND is sensitive to the z/g flux, 

^'e + P — ^ + e"*" , (6.1) 

from 17 reactors that are located reasonably close to the detector. The distances 
from the different reactors to the experimental site vary from slightly more than 
80 km to over 800 km, while the majority (roughly 80%) of the neutrinos travel from 
140 km to 215 km (see, e.g., for details). 

KamLAND "sees" the antineutrinos by detecting the total energy deposited by 
recoil positrons, which are produced via reaction (|6.1| ). The total visible energy, 

-^visible; is 

-^visible = Ee+ + me = Ep + TTlp - TJlr, + ITt^, (6.2) 

where Ee+ is the total energy of the positron and mp, mn,and rrie are, respectively, the 
proton, neutron, and electron mass. The positron energy is related to the incoming 
antineutrino energy by E(.+ = — 1.293 MeV, up to small corrections related to 
the recoil momentum of the daughter neutron. KamLAND is expected to measure 
the visible energy with a resolution better than a[E)/E = 10%/y/E, for E in MeV 

MM- 



In section ^]T] we present the predictions of the LMA solution (a) (cf. figure |l]a) 



for the charged current capture rate (eq. |6.1|) and for the distortion of the visible 
energy spectrum. We characterize the spectral distortion by the first and second 
moments of the energy distribution. No observable effect is predicted for the other 
currently allowed oscillation solutions. We describe in section |6]^ the calculational 
details of how the predictions were obtained. 

6.1 Predictions of global solution for charged current rate and spectrum 
distortion 



Table |g and Figure ^ summarize our principal predictions for the KamLAND reactor 
experiment. 

The reduced charged current ratio, [CC], given in Table ^ is the ratio of the 



observed event rate for eq. ( pTTP divided by the event rate predicted by the standard 



solar model (BPOO) and the standard electroweak theory (no oscillations). Thus 

[CC] = (6.3) 

If the predictions of the standard solar model and the standard electroweak model 
are both correct, then [CC] = 1.0. 
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Table 8: KamLAND event ratio and first and second visible energy moments. Here 
[CC] is the reduced charged current event ratio (defined below) and (5£'visibio and 5a are 
the first two moments of the visible energy spectrum. The results are tabulated for two 
different thresholds of the visible energy, 1.22 MeV (columns two through four) and 2.72 
MeV (columns five through seven) . The ranges listed correspond to the 1 a and 3 a allowed 
regions for our standard LMA oscillation solution, (a). 
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Figure 15: KamLAND reactor predictions. The figure shows the 3c7 allowed predic- 
tions for the CC and the first and second moments of the visible energy spectrum relative 
to the expectations based upon the standard solar model and the standard electroweak 
model (cf. Table |8|). The results are displayed for a visible energy threshold of 1.22 MeV 
(cf. eqs. and 6.5). 



We find that the allowed range of the charged current ratio computed for a 1.2 
MeV total visible energy threshold and the LMA region shown in figure |I|a is 



[CC] = QAAtr^ila) [0.44^^:f?(3a)], ^threshold = 1.22 MeV. (6.4) 

For a somewhat higher visible energy threshold in which the background is less 

I, Q), we find 



problematic (cf. refs. 

[CC] 



0.39l°:g^(la) [0.39+°i^(3a)], 



threshold 



2.72 MeV. 



(6.5) 



Table |^ and Figure represent the distortion of the visible energy spectrum by 
the fractional deviation from the undistorted spectrum of the first two moments of 
the energy spectrum (cf. the discussion in section |4.4| ) In the absence of oscillations, 
one expects (^vis)o = 3.97 MeV and ((t)o = 1.26 MeV for ^threshold = 1-22 MeV 
{{E^is)o = 4.33 MeV and {cr)o = 1.06 MeV for ^threshold = 1-72 MeV). The predicted 
distortion of the energy spectrum, which can be as large as 20%, may be more difficult 
to measure than the predicted deviation from unity of the reduced charge current 
rate, [CC]. If a significant distortion of the spectrum is observed, the magnitude of 
the distortion will severely restrict the possible range of the oscillation parameters. 

The results given here are in general agreement with other calculations that were 



made using previously determined allowed regions [11, B3, 54, 55, Bq, 57 



However, our calculations are the first that we know of that present the predicted 
distortions of the visible energy spectrum in terms of the first and second moments 
and which give the predicted values of the reduced charge current ratio for two 
separate energy thresholds (and for Icr and 3a deviations). 

6.2 Calculational procedures and details 

The antineutrino spectrum which is to be measured at KamLAND depends on the 
power output and fuel composition of each reactor (both change slightly as a function 
of time) and on the cross section for reaction ( |6.1|) . For the results presented here, 
we follow the fiux and the cross section calculations and the statistical procedure 
described in ref. 



55 



We define one "KamLAND-year" as the amount of time it 
takes KamLAND to observe 800 events with visible energy above 1.22 MeV. This is 
roughly what is expected after one year of running (assuming a fiducial volume of 1 
kton), if all reactors run at (constant) 78% of their maximal power output [^. We 
assume a constant chemical composition for the fuel of all reactors (explicitly, 53.8% 
of 235u, 32.8% of 239pu, 7.8% of 238U, and 5.6% of ^^^Pu, see |3). 



The shape of energy spectrum of the incoming neutrinos can be derived from a 



phenomenological parametrization, obtained in ^1 



(X e 



ao+aiE„+a2Ef, 



where the coefficients depend on the parent nucleus 
different isotopes we used are tabulated in [pT 



(6.6) 

The values of for the 
These expressions are good 



approximations of the (measured) reactor fiux for values of E^ >2 MeV. 
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The cross section for z/g + p — >■ e"*" + n has been computed including corrections 
related to the recoil momentum of the daughter neutron in ||6^. We used the hy- 
drogen/carbon ratio, r=1.87, from the proposed chemical mixture (isoparaffin and 
pseudocumene) |5^. It should be noted that the energy spectrum of antineutri- 
nos produced at nuclear reactors has been measured accurately in previous reactor 
neutrino experiments (see |^ for references). For this reason, we will assume for 
simplicity that the standard (without oscillations) antineutrino energy spectrum is 
known precisely. Some of the effects of uncertainties in the incoming flux on the 
determination of oscillation parameters have been studied in |^ and are estimated 
to be small. 

The calculation is first done for visible energies 1.22 < E^s < 7.22 MeV and in 
the computation of the shifts on the visible energy moments we assume 12 energy 
bins (binwidth is 0.5 MeV). There still remains, however, the possibility of irreducible 
backgrounds from geological neutrinos in the lower energy bins {Ems ^ 2.6 MeV) 
[0, ^ To verify how this possible background may affect the results we have 
repeated the analysis discarding the three lower energy bins i.e. considering only 
events with visible energies 2.72 < E^s < 7.22 MeV. 



7. Discussion and summary 

In this section, we summarize and discuss our principal results that are presented in 
the previous sections. 

We summarize in section ^[T] our conclusions regarding solar model predictions 



and present in section the currently allowed solar neutrino oscillation solutions 
that are found with three different analysis strategies. We describe the predictions of 
these allowed solutions for the SNO experiment in section [77^ and for the BOREXINO 



and KamLAND experiments in section 7.4. In section 7.5 we discuss what can be 



learned from the BOREXINO and KamLAND solar neutrino experiments, despite 
the fact that the favored LMA, LOW, and VAC experiments all predict similar event 
rates. We describe the unique sensitivity of the reactor KamLAND experiment in 
section [7.6| . Finally, we summarize in section |7.7| the robustness of the neutrino 
oscillation predictions and the possibilities for uniquely identifying the oscillation 
solution. 

We present predictions for both la and 3a ranges of the currently allowed neu- 
trino oscillation parameters. 

We use the notation [Q] to represent the measured or predicted value of the 
quantity Q divided by the value expected for Q if one assumes the correctness of the 
standard solar model (BPOO) and also assumes that nothing happens to the neutrinos 
after they are created in the sun. 
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7.1 Uncertainties in model predictions of solar neutrino fluxes 

How does the *B neutrino flux predicted by the solar model agree with the value in- 
ferred from the measurements by SNO for the CC flux and by Super-Kamiokande 
for the neutrino-electron scattering rate? Combining the SNO and Super-Kamiokande 
measurements, the SNO collaboration infers |0](see also ^) a total active 
neutrino flux (i/g + t'^i + ^r) of 0(^B; measured) = (5.44 ± 0.99)10^ cm~^s~^. Com- 
paring the measured ^B neutrino flux with the flux [|12[ calculated using the 1998 
standard value for ^17(0), 0(SB; BPOO) = 5.05(1^|]:?^)106 cm-^s'^, one finds 



0(SNO + Super — Kamiokande) 

r = 1.08 ± 0.20 (measurement) ± 0.18(solar model). 

0(BPOO) 

(7.1) 

Thus the measured flux agrees with the BPOO solar model flux to within O.So", com- 
bined experimental and theoretical errors. If we use instead the flux (see table j!]) ob- 
tained with the more recent (Junghans et al.) value for 5*17(0), 5.93 x 10~^ (l.OO^Qjg), 
we find 



0(SNO -|- Super — Kamiokande) 

-— — -— = 0.91 ± 0.17( measurement) ± 0.14(solar model). 

0(BPOO + New^B) ^ ' ^ ' 

(7.2) 

Therefore the measured flux agrees with the BPOO -|- New^B solar model flux to 
within 0.4cr, combined experimental and theoretical errors. In both cases, the mea- 
sured and the solar model fluxes agree to better than la. The agreement between 
the predicted and measured ^B solar neutrino flux is a significant confirmation of 
the standard solar model since the rare *B flux depends upon the 25th power of the 



central solar temperature [34 



We present calculated results in this paper for two values of the crucial low 
energy cross section factor, S'i7(0), namely, 5*17(0) = 19^.2^^^ [0] (the 1998 standard 
value) and 6*17(0) = (22.3 ± 0.9)eVb (the most recent and precise measurement). 
These two values of 6*17(0) represent well the range of measured cross section factors 



obtained in recent experiments [T^. The reader can therefore see at a glance the 
dependence of the calculated quantities upon the adopted value of the cross section 
^17(0). 

Essentially all global analyses of solar neutrino oscillations make use of neutrino 
fluxes and their uncertainties obtained from a solar model. If the eight most im- 
portant solar neutrino fluxes(cf. table |l| ) were allowed to vary without constraints, 
there would be too many free parameters to make possible an efficient study of solar 
neutrino oscillations. If the BPOO model is adopted in the analysis, then the total 
p — p flux is assumed known to ±1% and the total "^Be flux is assumed known to 
±10%. 
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Guidance from a solar model is required just to decide what are the most im- 
portant neutrino fluxes to include in the analysis. For example, all current analyses 
of solar neutrino oscillations include a subset of the eight fluxes calculated in the 
standard solar model and listed in Table |l] but neglect the fluxes from the reactions 
e- + 3He -> ^H + z/e, + ^ ^Be + u^, e- + '^^C + u^, e" O ^ ^^N + z/^, 

and e~ + ^''F — » ^^O + z/g. One must use the parameters of a solar model to show that 
the neglected fluxes are negligible and that the included fluxes are important. 

Even for a Bayesian analysis (see ref. ||66|), one needs guidance from a solar 
model to determine which fluxes should be included and which should be neglected 
and, e. g., to decide if the priors should be expressed in terms of lO^cm^s"^ for the 
^Be flux and 10® cm^s~^ for the ®B flux or vice versa. Nevertheless, the data for solar 
neutrino experiments is becoming sufficiently precise and constraining that Bayesian 
analyses can provide important and complimentary insights to results obtained with 
other statistical techniques (see, e.g., figure 20 of ref. p6|] ). 

In a few cases, authors have carried out analyses in which the *B neutrino flux, 
but none of the other fluxes, is completely unconstrained by the solar model pre- 
dictions, which is our case (c) in section |^ and in figure |l|c. The experimental data 
from solar neutrino experiments are not yet sufficient to permit a very restrictive 
analysis for unknown neutrino oscillation parameters if neutrino fluxes other than 
the ®B neutrino flux are allowed to be free variables. 

The undesirable but currently unavoidable practice of using solar models to help 
constrain neutrino parameters is, ironically, the opposite practice to what was envi- 
sioned in the early days of solar neutrino research. The original chlorine experiment 
was proposed to test solar models using the assumed well-known properties of neu- 
trinos [|67| . 

In order to document what we use in this paper and to facilitate analyses of 
solar neutrino data by other groups, we summarize in section H the uncertainties 
and the best-estimates for solar neutrino fluxes derived from the BPOO solar model. 
We describe the results that are obtained when adopting the precise Junghans et al. 
value of S'i7(0) and also the results that are obtained using the 1998 standard value 
of 5*17(0) (see especially table |l| and table || and related comments in the text). 

We stress the continued importance of precise measurements of the low energy 
cross section factor, S'i7(0), for the '^Be(p,7)*B reaction. The neutrino fluxes mea- 
sured in the Super-Kamiokande, SNO, and ICARUS solar neutrino experiments all 
depend linearly upon this cross section factor and the standard model prediction of 
the event rate in the chlorine experiment is also dominated by the ^B neutrino flux. 

The uncertainty in the low energy cross section factor, S'34(0), for the ^He(^He,7) 
^Be reaction dominates the uncertainty in the solar model calculation of the ^Be solar 
neutrino flux. The total uncertainty in the solar model calculation of the ^Be neutrino 
flux is 9% and the ^He(^He,7)^Be reaction contributes 8% to the total uncertainty 
that is computed by quadratically combining uncertainties from different sources. 



- 39 - 



The cross section factor 5*34(0) is also the largest nuclear physics uncertainty in the 
prediction of the solar neutrino flux if one adopts the Junghans et al. value for 
the uncertainty in 5*17(0). 

Now that there are precise measurements of 5i7(0) underway and completed, we 
believe the highest priority nuclear astrophysics measurement for the future is the 
precision determination of the low energy cross section factor for the ^He(^He,7)^Be 
reaction. A measurement of 5*34(0) to an accuracy of better than 4% is necessary 
in order to decrease the uncertainty in the reaction ^He(^He,7) ''Be to where it is no 
longer the largest uncertainty in the prediction of the ^Be solar neutrino flux. 

7.2 Global neutrino oscillation solutions 

Using three different analysis strategies that span the range of previously used strate- 
gies, we determine the globally allowed solar neutrino oscillation solutions that are 
consistent with all the available solar and reactor data. The results are summarized 
in figure the calculations on which this figure is based used the new Junghans et 
al. value of 5i7(0). Table ^ gives the best-fit values of Am^ and tan^ 6*, as well as 
the local value of Xmm f*^^ each oscillation solution; the results presented in table ^ 
were obtained using our standard analysis strategy in which we take account of the 
BPOO -|- New ^B predicted fluxes and estimated uncertainties and include the Super- 
Kamiokande day and night energy spectrum but not the total rate. The 3a upper 
limit to the allowed value of Am^ lies between 3.0 and 7.5, in units of 10""^ eV^, 
depending upon what is assumed about 5i7(0) and the analysis strategy. 

The LMA solution is favored, but the LOW and VAC solutions are also allowed at 
a CL. corresponding to 3a. Other solutions such as oscillations into sterile neutrinos 
or the SMA solution for active neutrinos are disfavored at 3a if the standard analysis 
strategy is adopted, but SMA is barely allowed at 3a if the ^B neutrino flux is 
unconstrained by the solar model predictions and uncertainties. 

Figure |] is a "Before and After" comparison of the effect of the ^B production 
cross section on the global oscillation solutions. The only difference between the 
"Before and After' panels in figure ^ is the replacement in the analysis of the 1998 
standard value of 5*17(0) = 19t2eVb ||19[ by the value obtained by the recent precise 
measurement, i.e., 5i7(0) = (22.3 ± 0.9) eVb §. Thus the "before" (left) panel of 
figure 1^ corresponds to the results shown in the left panel of figure 9 of ref. |T^. The 
change in the value of 5*17(0) is sufficient to drive the SMA and Just So^ solutions 
over the edge of the allowed region; they are allowed at 3a with the older value of 
5*17(0) but not with the newer value. 

The global oscillation analysis yields a 3a allowed range for the inferred total 
^B solar neutrino flux expressed in terms of the best-estimate flux predicted by the 
BPOO model with the Junghans et al. value of 5*17(0). From table H, we find for 
active neutrinos 

0.40 < /b < 1.36, (7.3) 
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for the case in which the neutrino flux is unconstrained in the analysis. The 
best-fit value of /b = 1.0 for this unconstrained case. The 3a allowed range found 
here is slightly smaller than determined directly for active neutrinos by the SNO 
collaboration by comparing the SNO and Super-Kamiokande results. The SNO col- 



laboration found 10] /b = 0.92 ± 0.50, 3a range. We show in section ^]2| that if one 
considers an arbitrary admixture of sterile and active neutrinos, the 3a upper limit 
of /b increases to 2.9. 

All of the global analyses of solar neutrino experiments that include the im- 
portant Super-Kamiokande data on the electron recoil energy spectrum use the 
many energy bins provided by the Super-Kamiokande group (19 energy bins in the 
last report). It would be very instructive to be able to carry out global analyses 
while representing the spectrum by only one or two measured quantities, the first or 
the first and second moments. One could then determine the change in the global 
solutions that result from giving the measurement of the energy spectrum the same 
prominence in the analysis as one or two measurements of the total rates^. 

We describe in Section the technical differences between the three analysis 
strategies. This section is intended primarily for neutrino analysis junkies. 

Table |] presents the contributions of individual sources and the total rates pre- 
dicted by the favored neutrino oscillation schemes, LMA, LOW, and VAC, for the 
chlorine and the gallium radiochemical experiments. The predictions are given in 
the table for our best-fit solutions obtained using the standard analysis strategy (a). 
The two measured rates are also listed. The errors are dominated by the solar model 
uncertainties for analysis strategy (a), which is apparent from table ^by comparing 
the calculated and measured values. The theoretical uncertainties expressed in SNU 
for the chlorine and gallium rates are greatly reduced for the oscillation solutions 
(table ^ compared to the no-oscillation values (table [1|). However, the theoreti- 
cal uncertainties expressed as percentages of the total rates are comparable for the 
oscillation and no-oscillation scenarios. 

All three of the oscillation scenarios yield total event rates in good agreement 
with the measured values for the gallium experiments. However, the calculated rates 
for the LOW and the VAC solutions are in poor agreement with the measured values 
(discrepancies of 2.7 and 3.8 a, respectively). The global VAC solution is acceptable 
only because one can choose the normalization such that the predicted distortion 
is small for the recoil energy spectrum measured by Super-Kamiokande. Of course. 



^We already know from previously published calculations (see refs. [|l3, 15, 2T| that are significant 
differences in the allowed regions between the extreme cases in which only the total rates of the 
solar neutrino experiments are considered and the case in which all 19 of the Super-Kamiokande 
energy bins are included in addition to the total rates. In order to properly take account of 
the characteristics of the detector, which can influence the moments 37, 0, the experimental 



collaboration must publish both the measured moments and the values expected for the moments 
if there is no distortion of the spectrum. 
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Source 


CI 


Ga 


CI 


Ga 


CI 


Ga 




(SNU) 


(SNU) 


(SNU) 


(SNU) 


(SNU) 


(SNU) 




LMA 


LMA 


LOW 


LOW 


VAC 


VAC 


DP 

irsr 





41.8 





38.7 





44.3 




0.12 


1.49 


0.11 


1.35 


0.16 


1.95 




n ni 






u.uo 


u.uo 


n OA 


^Be 


0.62 


18.7 


0.58 


17.5 


0.54 


16.4 


«B 


2.05 


4.27 


2.94 


6.13 


3.95 


8.33 


13N 


0.05 


1.80 


0.04 


1.69 


0.05 


2.01 


150 


0.17 


2.77 


0.16 


2.64 


0.20 


3.29 


17p 


0.00 


0.03 


0.00 


0.03 


0.00 


0.04 



Total 3.01 ±0.31 71.0 ±2.6 3.85 ± 0.45 68.1 ±2.8 4.93 ± 0.60 76.4 ± 3.2 

Table 9: Neutrino oscillation predictions for the chlorine and gallium radiochemical ex- 
periments. The predictions are based upon the global analysis strategy (a), described in 
section ^ and make use of the neutrino fluxes given in table |l] for the BPOO solar model, 
the neutrino absorption cross sections given in ref. and the Junghans et al. value 
of 5'i7(0) The rates and la uncertainties are presented for the best-fit values of the 
allowed solutions listed in table |^. The total rates should be compared with the standard 
solar model values of table |l|, which are, 8.65^J'^ (chlorine) and 130^y (gallium), and the 
measured values, 2.56 ±0.23 (chlorine, see ref. ||l|)and 75.6 ± 4.8(gallium, see ref. ||2|, |3|, Q). 



for strategy (c), in which the ^B neutrino flux is a free parameter, one can obtain 
global allowed solutions that are much better fits to the chlorine rate. 

The predicted rates for the Super-Kamiokande experiment (in units of x 10^ 
cm-2s-i) are: 2.39l[]:| (LMA), 3.02toil (LOW), and 3.8i;[J|?(VAC), which should 
be compared with the measured rate of (2.32 ± 0. 03^9 07) x 10^cm~^s~^ Q. For the 
SNO CC measurement, the predicted rates are (in units of xl0^cm~^s^^) : 1.72^q;26 
(LMA), 2.50l°:| (LOW), and 3.261°;^^ (VAC), which should be compared with the 
measured rate of (1.75 ± 0.07+°;^^) x lO'^cm-^g-i 0. 

7.3 Predictions for SNO 

All three analysis strategies yield essentially the same 3a range for the neutral current 
to charged current double ratio(see figure ||) predicted by the favored LMA, LOW, 
and VAC solutions, namely, 1.4 < [NC]/[CC] < 5.3(For analysis strategy (c), with 
uses an unconstrained '^B neutrino flux, the upper limit extends to 6.2.) The la 
predicted range is 3.5 ± 0.6 for our standard analysis strategy. 

The predicted range of [NC]/[CC] for the favored oscillation solutions is also 
rather insensitive (see figure ^) to the choice of the ^B production cross section within 
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the range of the 1998 standard value, ref. [19|, and the recent precision determination, 
ref. P]. If the is allowed to vary freely, our analysis strategy (c), then the SMA 
solution for active neutrinos is allowed at the 3a CL. and for this solution 1.15 < 
[NC]/[CC] < 1.31. Figure | shows that all of the values of [NC]/[CC] predicted by 
the allowed neutrino oscillation solutions are separated from the no oscillation value 
of [NC]/[CC] = 1.0 by more than the expected experimental uncertainty. 

For the average difference in the CC day-night effect, ylN_D(SNO CC), all three 
analysis strategies also yield essentially the same results, as is shown in figure ^ 
Figure H shows that the predictions for y4N_D(SN0 CC) are also robust with respect 
to likely changes in the value of 5*17(0). The 3a range in percent for all the oscillation 
solutions is 0.0 < An-d(SNO CC) < 0.21 (The la range is An-d(SNO CC) = 
8.15 ±5.15 % .) 

The predicted average difference in the day-night effect for neutrino-electron 
scattering, An-d(SNO ES), depends upon which among the currently allowed global 
oscillation solutions (LMA, LOW, VAC, or SMA) is assumed (see figure 0). All 
of the results are bounded by the 3a allowed range for Super-Kamiokande, i. e., 
An-d(SNO ES) lies between 0% and 11%. Figure || shows the predicted correlation 
between the day-night effect for the CC and the day-night effect for neutrino-electron 
scattering. The predicted correlation between 74n-d(SNO ES) and y4N_D(SN0 CC) 
constitutes an important consistency check for the oscillation solution. 

The first and second moments of the electron recoil energy spectrum from CC 
reactions summarize most of the measurable information about the energy spectrum 
of neutrinos that are observed by SNO (and by Super-Kamiokande). The low 
order moments may represent a more appropriate way to characterize the relatively 
undistorted energy spectrum rather than to provide event rates in a large number 
of separate spectral energy bins(cf. ref 0). SNO has a significant advantage over 
Super-Kamiokande in measuring potential distortions of the energy spectrum due to 
neutrino oscillations, because the electron is the only light particle (hence it takes 
most of the energy) in the CC reaction, eq. ( |4.1| ). 

The current set of allowed oscillation solutions do not predict distortions that 
are large enough to be detected by SNO at a high level of significance. Figure | and 
figure |TU|, when combined with the discussion in section |4.4] of the likely experimen- 
tal uncertainties, establish this pessimistic prediction. Nevertheless, it is important 
to measure accurately the first two moments of the recoil energy distribution. The 
pessimistic prediction summarized here could be wrong if the correct oscillation so- 
lution is one not favored by the initial pioneering set of solar neutrino experiments 
(see further discussion of this point in section fLTl ). 

If we are lucky, if Nature is kind, then one measurement may define rather 
well the solar neutrino oscillation parameters. Figure |ll] shows that only relatively 
small regions in neutrino oscillation parameter space predict that SNO will measure 
either [NC]/[CC] > 3.3 or An-d(SNO CC) > 0.1, If either of these inequalities is 



- 43 - 



established experimentally, then the oscillation parameters Am^ and tan^ 9 will be 
strongly constrained. 



7.4 Solar neutrino predictions for BOREXINO and KamLAND 

Table |^ and figure |T3| show that the currently LMA, LOW, and VAC oscillation solu- 
tions predict similar values for the neutrino-electron scattering rate for ^Be neutrinos, 
pBe] = 0.655 ±0.035, 1(t(0.66 ± 0.13, 3a). The SMA solution, which is allowed at 3a 
only in the "free analysis strategy (cf. figure |l|c), predicts a much smaller value, 
f Be] < 0.34 . 

Figure |l^ and table |^ show that the LOW solution is the only currently allowed 



(at 3a) neutrino oscillation solution that predicts a significant day-night variability 
in BOREXINO or KamLAND. If a difference between day and night rates is actually 
observed, this will very strong evidence in favor of the LOW solution. 

We have tested the robustness of the predictions for the BOREXINO and Kam- 
LAND experiment in several ways: i) we computed the predicted results for all three 
of the analysis strategies described in section ii) we compared the results at the 
terrestrial locations of BOREXINO and KamLAND; and iii) we compared the pre- 
dictions with and without including events from neutrino sources other than ^Be. In 
all cases, the variations are small (and are given quantitatively in section §). 

How well can we predict [^Be] if SNO measures either [NC]/[CC] > 3.3 or 
y4N-D(SN0 CC) > 0.1? Unfortunately, [''Be] is not a unique function of the neu- 
trino mixing parameters, Am^ and tan^ 9. Instead, [''Be] can take on the same value 
for many different pairs of {Am?, iai?9). The predicted range of [^Be] over the full 
LMA region allowed at 3a is from 0.585 to 0.76. For the parameters corresponding to 
[NC]/[CC] > 3.3, the range is almost as large: 0.585 to 0.74. For An-d(CC) > 0.1, 
the predicted range of [^Be] is 0.585 to 0.72. Thus a hypothetical (and optimistic) 
experimental determination of the behavior of the survival probability, represented 
by figure 0, in the 6 MeV to 10 MeV region of neutrino energies most effectively 
probed by SNO would not be sufficient to allow a precise prediction of the survival 
probability at 0.86 MeV, which corresponds to the energy of the ^Be neutrino line. 

7.5 What will we learn from BOREXINO and KamLAND solar neutrino 
experiments? 

At first impression, it may seem somewhat discouraging that the currently most 
favored oscillation solutions, LMA, LOW, and VAC, all predict very similar event 



rates for neutrino-electron scattering in BOREXINO and KamLAND (see figure [T^). 
However, this is really an advantage. 

If the currently favored oscillation solutions and the standard solar model are 
correct, then BOREXINO or KamLAND must confirm that [^Be] = 0.66 ± 0.13. 
Thus the fact that the predicted event rates for BOREXINO and KamLAND are 
similar for the LMA, LOW, and VAC solutions means that a measurement of the 
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neutrino electron scattering rate is a critical test of the validity of the standard 
solar model prediction and the favored oscillation solutions. In addition, we may be 
somewhat surprised and the measurements may show that [''Be] < 0.34, favoring an 
SMA solution. 

Also, a large measured value for the day-night effect in BOREXINO or Kam- 
LAND would imply the correctness of the LOW oscillation solution, which is cur- 
rently allowed but is less favored than the LMA solution (cf. Table 

Experiments like BOREXINO and KamLAND that detect neutrinos with ener- 
gies less than 1 MeV are necessary to test the validity of solar neutrino oscillation 
solutions. The solutions explored in this paper, and in related papers by other au- 
thors, are primarily constrained by data that refer to the relatively high energies(> 5 
MeV) for solar neutrinos to which the chlorine, Kamiokande, Super-Kamiokande and 
SNO experiments are primarily or exclusively sensitive. Of all the solar neutrino ex- 
periments performed so far, only the gallium experiments, SAGE and GALLEX -|- 
GNO, have a large sensitivity to low energy neutrinos. The gallium experiments are 
radiochemical experiments and therefore do not measure neutrino energies. 

BOREXINO and KamLAND will also tell us something new and critical about 
the Sun. We have very little direct observational information about the important 
^Be solar neutrino flux. If one supposes that ^Be is the only source contributing to 
the observed rates in the chlorine and gallium experiments, then the the Sex upper 
limit implied by the chlorine measurements [|l| is 2.8 times the BPOO predicted flux. 
The upper limit implied by the gallium experiments 0, 0, ^ is 2.6 times the BPOO 
predicted flux. All of the existing experiments are consistent with a ^Be flux that is 
identically zero. Thus the direct experimental constraints only require that 

0.0 < l\ ' < 2.6, 3a. 7.4 

0(7Be)BPoo 

It is important to test experimentally whether the solar model prediction for the ''Be 
neutrino flux is correct. The branch of the p — p fusion chain that is represented by 
^Be neutrinos occurs in 15% of the completions of the chain, according to the BPOO 
solar model. 

In order to test the solar model prediction of the ^Be total neutrino flux, we must 
perform at least one additional experiment beyond the BOREXINO and KamLAND 
u — e scattering experiments, which measure a linear combination of CC and NC 
events. Either a CC measurement or a NC measurement of the ''Be line is essential 
to test the solar model in a way that is completely free of all influence of solar models. 

7.6 Predictions for the KamLAND reactor experiment 

The KamLAND reactor experiment may provide definitive evidence for or against the 
LMA oscillation solution. The reduced charged current rate is predicted by the LMA 
solution to satisfy [CC] = O.Utof^ila) [0.44+°:29(3c7)], ^threshold = 1-22 MeV, with 
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a somewhat larger allowed range for a higher energy threshold. Table ^ summarizes 
the predicted distortion of the energy spectrum (in terms of the first and second 
moment) as well as the rate of the CC reaction for two plausible energy thresholds. 
If the LMA solution is correct, then the KamLAND experiment should observe a 
measurable deficit in the charged current rate. The spectral distortion, if measured, 
will determine the solar neutrino oscillation parameters with unprecedented precision. 

7.7 Robustness and Uniqueness 

The neutrino oscillation solutions that describe all the available solar and reactor data 
are, with one exception, robust with respect to variations among the three analysis 
strategies we have used in this paper. The SMA solution is the sole exception, and 
only just barely an exception. This solution exists at the 3a CL. if the neutrino 
flux is treated as a free parameter (cf. figure but the SMA is absent for all three 
analysis strategies at 99% C.L. 

The predictions for the quantities that are expected to be most easily measured 
by the SNO, BOREXINO, and KamLAND experiments are also robust with respect 
to analysis strategies. 

If we are lucky, one of these experiments may determine uniquely the solar 
neutrino oscillation parameters. SNO could potentially identify the LMA solution as 
correct by observing a large value for either the neutral to charged current ratio or the 
day-night difference. KamLAND could also establish the LMA solution by finding 
an appreciable deficit of reactor antineutrinos in the first phase of the experiment. 
On the other hand, BOREXINO and KamLAND could definitively select the LOW 
solution if a large day-night effect is observed. 

All of the above comments depend upon the validity of the global theoretical 
analysis which fits oscillation solutions to the reported results of the first six so- 
lar neutrino experiments: chlorine, Kamiokande, SAGE, GALLEX + GNO, Super- 
Kamiokande, and SNO. Because the analyses depend upon experimentally untested 
constraints implied by the standard solar model and because there is insufficient re- 
dundancy in the experiments performed so far, there may be startling surprises when 
solar neutrino oscillations are probed in the future in different experimental ways. 

7.8 Late breaking news from Super-Kamiokande 

The Super-Kamiokande collaboration has discussed preliminary data for 1496 
days of observation, a 19% increase in the length of data-taking over the previously 
reported data set of 1258 days of observation. We are grateful to M. Smy for making 
the preliminary data available to us so that we could access the robustness of the 
predictions made in this paper to a modest amount of additional data. 

Figure ^ shows the global solutions for the same three analysis strategies and 
input data as were used in producing figure 0, except that figure ^ makes use of 
1496 days of Super-Kamiokande data. At first glance, figure |T] and figure |TB| are very 
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Figure 16: Global oscillation solutions including 1496 days of Super-Kamiokande 
observations. Figure 16 is the same as figure |l] except that data from 1496 days of 
Super-Kamiokande observations were used in constructing figure |l^ whereas 1258 days of 
Super-Kamiokande observations were used in constructing figure |l[ 



similar to each other. With the new data, there is a small reduction in the lower-mass 
region of the LMA solution and in the upper-mass region of the LOW solution. The 
primary reasons for the changes that do occur are the somewhat smaller differences 
between the day and the night data and the slightly flatter (with respect to the 
undistorted standard spectrum) recoil energy spectrum. We spare the reader the 
details of the best-fit points (which are shown in figure |T6[) since these values bounce 
around within the allowed region as new data are obtained. 

At what CL is LMA the unique solution to the solar neutrino data? In the 
global solution (strategy a) represented in figure |1|, which includes 1258 days of 
Super-Kamiokande data, LMA is the only solution at 96% {2a) CL and in figure |TB|, 
for 1496 days of Super-Kamiokande observations, LMA is the only solution at 97% 
(2.1(t) CL. Even before the availability of the SNO data, LMA was the only solution 
at 86% (1.5(t) CL [|^. These results are in good agreement with those described 
by Smy ||68| (who analyzed the Super-Kamiokande data with more zenith-angle bins 
and fewer spectral energy bins than we do). We conclude that LMA is the preferred 
solution with a CL that has been increasing slowly with time. 



Table |T0] presents the la and 3a ranges for the principal solar neutrino observ- 
ables that were calculated in sections In constructing table |10|, we again use 



analysis strategy (a) and the same input data as in the previous calculations ex- 
cept that we now include 1496 days of Super-Kamiokande data.It is instructive to 
compare the entries in table [1^ with the previously tabulated results. The changes 
due to the new data are smaller than or comparable to the differences among the 
predictions of the three analysis strategies for the total event rates and the energy 
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Table 10: Predictions including 1496 days of Super-Kamiokande data. This table 
presents the la and 3cr ranges for the measurable quantities calculated in sections |^]6|. 
We use here analysis strategy (a) and the same input data as in the previous calculations, 
except that 1496 days of Super-Kamiokande data were included. Results are presented for 
measurables in the SNO, BOREXINO, and KamLAND experiments. The threshold of the 
recoil electron kinetic energy used in computing the SNO observables for this table is 6.75 
MeV and for the KamLAND reactor observables, £'th = 1.22 MeV. The recoil energy range 
for the BOREXINO experiment is the same as adopted in section ^. 

shifts, i. e., for the quantities [NC]/[CC], 5T(SN0 CC), fBe], [CC] (KamLAND), 
and 6E^is (KamLAND). We conclude, once again, that these predictions are robust. 

The predicted maximum value at 3cr of the day-night asymmetry is reduced by 
about 15% for the CC measurement by SNO (from 21% to 18%) and by about 10% 
for BOREXINO (from 41% to 38%). These changes reflect the difficulty in measuring 
day-night asymmetries, which are determined by relatively small differences between 
two large numbers. 
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